FUNCTIONAL CALCULUS OF DIRAC OPERATORS AND 
COMPLEX PERTURBATIONS OF NEUMANN AND DIRICHLET 
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Abstract. We prove that the Neumann, Dnichlet and regularity problems for 
divergence form elliptic equations in the half space are well posed in L2 for small 
complex Loo perturbations of a coefficient matrix which is either real symmetric, 
of block form or constant. All matrices are assumed to be independent of the 
transversal coordinate. We solve the Neumann, Dirichlet and regularity problems 
through a new boundary operator method which makes use of operators in the 
functional calculus of an underlaying first order Dirac type operator. We estab- 
lish quadratic estimates for this Dirac operator, which implies that the associated 
Hardy projection operators are bounded and depend continuously on the coeffi- 
cient matrix. We also prove that certain transmission problems for fc-forms are 
well posed for small perturbations of block matrices. 
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1. Introduction 

In this paper we prove that the Neumann, Dirichlet and regularity problems are 
well posed in L2(R") for divergence form second order elliptic equations 

(1.1) diYt,,A{x)Vt,,U{t,x) = 

on the half space R""*"^ := {{t,x) G R x R*^ ; t > 0}, n > 1, when A is a small 
complex Loo perturbation of either a block matrix, a constant matrix or a real sym- 
metric matrix. Furthermore, the matrix A = (aij(a;))"j=o ^ -^oo(R"'; '^(C'^"'"^)) is 
assumed to be t-independent with complex coefficients and accretive, with quanti- 
tative bounds ll^lloo and ka, where /tyi > is the largest constant such that 

Re{A{x)v,v) > ka\v\\ for all v G C"+\ x G R". 

We shall approach the equation (11. ip from a first order point of view, rewriting it 
as the first order system 

' diYt,xA{x)F{t,x) =0, 
cm\t^^F{t,x) = 0, 



1.2) 



where F(t, x) = Vt,xU (t, x). Recall that a vector field F = FqCq + FiCi + . . . + F„e„ 
can be written in this way as a gradient if and only if curlj = 0, by which we 
understand that djFi = diFj, for all i, j = 0, . . . ,n. We write {eo, ei, . . . , Cn} for 
the standard basis for R""*"^ with eo upward pointing into R""*"^, and write t = xq 
for the vertical coordinate. For the vertical derivative, we write do = dt- Denote 
also by Fy := FiCi + . . . + F„e„, the tangential part of F, and write curlajFy = if 
djFi = diFj, for all i, j = 1, . . . ,n. 

In the formulation of the boundary value problems below, we assume that A = 
A{x) is a given coefficient matrix with properties as above. Furthermore, by saying 
that Ft{x) = F{t,x) satisfies we shall mean that Ft G Ci(R+; L2(R"; C^+i)) 
and that for each fixed t > 0, we have diVx{AF)^^ = —{A{x)doF)o, VxFq = OqF^^ and 
curlajFy = 0, where the derivatives on the left hand sides are taken in the sense of 
distributions. If this holds for F, then in particular we can write F = Vt,xU, with 
U G W2 ioc(^V~^) ^ ^^^^ ^ satisfies (11.11) in the sense that 



j I {A{x)Wt,xU{t, x), WtMt^ ^)) dtdx = 0, for all if G C^iYil^^). 

THE NEUMANN PROBLEM (Neu-A). 

Given a function 0(a;) G L2(R";C), find a vector field Ft{x) = F(t,x) in R"+^ 
such that Ft G C\R+; L2(R"; 0"+^)) and F satisfies for t > 0, and furthermore 
limj^oo Ft = and limt^o -^t = / in -^2 norm, where the conormal part of / satisfies 
the boundary condition 

n 

eo ■ (Af) = J2 Aojfj = 0, on R" = dRl+\ 

j=0 

For U, this means that the conormal derivative ^^(0,x) = 0(x) in L2(R"). 

THE REGULARITY PROBLEM (Reg- A). 

Given a function ip : R" — C with tangential gradient Vx'ip G L2(R"; C"), find 
a vector field Ft{x) = F{t,x) in R![+^ such that Ft G C\R+; L2(R"; 0"+^)) and F 
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satisfies f ll.2p for t > 0, and furtliermore linif^oo -^t = and limt-,oFt = / in L2 
norm, wliere tlie tangential part of / satisfies tlie boundary condition 

/„ = /lei + . . . + /„e„ = V.^, on R" = dRl+\ 

For U, tliis means tliat V^f/(x, 0) = V^?/'(a;), i.e. U{x, 0) = in W^{R''). 

THE DIRICHLET PROBLEM (Dir-A). 

Given a function u{x) G L2(R"; C), find a function f/t(a;) = U{t,x) in R""^""^ sucfi 
tliat Ut e C2(R+;L2(R";C)), Vt,^f/t G i^2(R"; C"+^)) and V*,^?/* satisfies 

(11. 2p for t > 0, and furtliermore lim^^oo Ut = 0, lim^^oo ^t,xUt = and limt_,o Ut = u 
in L2 norm. 

We shall also use first order methods based on (11. 2p to solve the Dirichlet problem. 
However, here we use a different relation between (II. ip and (II. 2p which we now de- 
scribe. Assume F{t,x) = J2k=o^i(^^-'^)^i ^ vector field satisfying (II. 2p . Applying 
dt to the first equation and using that curlf .^.F = yields 

= dt{diYt,,A{x)F) = diYt,.A{x){dtF) = dwt,.A{x){Vt,M, 

since the coefficients are assumed to be t-independent. Thus the normal component 
U := Fq satisfies (II. ip . Note that when A = I, the functions Fi, . . . ,Fn are conju- 
gates to U in the sense of Stein and Weiss [25]. From this we see that solvability 
of (Dir-A) is a direct consequence of solvability of the following auxiliary Neumann 
problem. 

THE NEUMANN PROBLEM (Neu^-A). 

Given a function (j){x) E L2(R'^; C), find a vector field Ft{x) = F{t,x) in R""^^ 
such that Ft G C\R+; L2(R"; 0"+^)) and F satisfies for t > 0, and furthermore 
lim^^oo Ft = and lim^^o Ft = f in L2 norm, where the normal part of / is Cq ■ / = 
/o = 0. 

The main result of this paper is the following Loo perturbation result for the 
boundary value problems. 

Theorem 1.1. Let Aq{x) = ((ao)ji(3;))"j=o t-independent, complex, accretive 
coefficient matrix function. Furthermore assume that Aq has one of the following 
extra properties. 

(b) Aq is a block matrix, i.e. (ao)oi(a^) = (a.o)jo(a^) = for all 1 < i < n and all 
X G R'^. 

(c) Aq is a constant coefficient matrix, i.e. Aq{x) = AQ{y) for all x, y E R*^. 
(s) Aq is a real symmetric matrix, i.e. (ao)ij(x) = (ao)ji(a;) G R for all < 

i,j<n and all x G R". 

Then there exists e > depending only on ||Ao||oo; the accretivity constant and 
the dimension n, such that if A E Loo(R"; /^(C"'"'"^)) is t-independent and satisfies 
\\A — v4o||oo < then the Neumann and regularity problems (Neu-A), (Neu^-A) 
and (Reg-A) above have a unique solution F{t,x) with the required properties for 
every boundary function g{x), being (f){x) and Vx'ip{x) respectively. Furthermore 
the Dirichlet problem (Dir-A) above has a unique solution U{t,x) with the required 
properties for every boundary function u{x). 
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The solutions depend continuously on the data with the following equivalences of 
norms. If we define the triple bar norm IWGtlW^ := WGtW^t^^dt and the non- 
tangential maximal function 

iV,(F)(x) := sup (-[ / \F{s,y)\Usdy) ' , 

where := \E\~^ and cq G (0,1), ci > are constants, then for the Neumann 
and regularity problems we have 

hh ~ II/II2 ^ sup WFth ~ WltdtFS ^ ||iV,(F)||2, 

and for the Dirichlet problem we have 

hh ~ sup WUth ^ III Willi ^ |||tV,f/i||| ^ \\N,{U)\\2. 

Moreover, the solution operators Sa, being SA^g) = F or Sa{u) = U respectively, 
depend Lipschitz continuously on A, i.e. there exists C < 00 such that 

\\Sa2 ~ Sai\\l2{'R^)^X < C\\A2 — ^i||loo(R.") 

when \\Ai — v4o||oo < ^ = 1;2, where ||-F||a' or \\U\\x denotes any of the norms 
above. 

Throughout this paper, we use the notation X ^ Y and X < Y to mean that 
there exists a constant C > so that X/C <Y< CX and X < CY, respectively. 
The value of C varies from one usage to the next, but then is always fixed. 

Let us now review the history of works on these boundary value problems, start- 
ing with the case of matrices of the form A which we now describe. By standard 
arguments. Theorem 11.11 also shows well posedness of the corresponding boundary 
value problems on the region Q above a Lipschitz graph E = {(t, x) ; t = g{x)}, 
where g : R" — R is a Lipschitz function. Indeed, if the function U {t, x) satisfies 
divt,xA{x)'Vt,xU = Omfl then U{t,x) := U{t + g{x),x) satisfies divt,xA{x)'Vt,xU = 
in Ri+\ where 



A{x) 



1 -{VMx)Y 
I 



A{x) 



1 
-V.(7(x) / 



Thus Theorem II . II gives conditions on A for which the Neumann problem (conormal 
derivative CQ-AVt^xU = {eo—'Vxg)-A'Vt,xU given), the regularity problem (tangential 
gradient VxU = {dtU)V xQ + V^f/ given), and the Dirichlet problem {U = U given) 
are well posed. Note that A is real symmetric if and only if A is, but that A 
being constant or of block form does not imply the same for A. For the Laplace 
equation A = J in f2, solvability of (Neu-J) and (Reg-/) was first proved by Jerison 
and Kenig (TU], and solvability of (Dir-J) was first proved by Dahlberg [12]. Later 
Verchota [2S] showed that these boundary value problems are solvable with the layer 
potential integral equation method. 

For general real symmetric matrices A, not being of the " Jacobian type" I above, 
the well posedness of (Dir-A) was first proved by Jerison and Kenig [17J, and (Neu- 
A) and (Reg- A) by Kenig and Pipher [21]. These results make use of the Rellich 
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estimate technique. For the Neumann and regularity problems, this integration by 
parts technique yields an equivalence 

(1-3) IISILaCR") ~ 

which is seen to be equivalent with the first estimate ||/|| ~ \\g\\ in the theorem 
above, and shows that the boundary trace / splits into two parts of comparable 
size. 

Turning to the unperturbed case where A = Aq satisfies (b), then (11. 3p is still 
valid, but the proof is far deeper than Rellich estimates. In fact, it is equivalent 
with the Kato square root estimate proved by Auscher, Hofmann, Lacey, Mcintosh 
and Tchamitchian in [3]. (For the non divergence form case ago 7^ 1, see pi|.) For 
details concerning this equivalence between the Kato problem and the boundedness 
and invertibility of the Dirichlet-to-Neumann map ^— ^ we refer to Kenig [T^ 
Remark 2.5.6], where also many further references in the field can be found. 

We now consider what is previously known in the case when A does not satisfy 
(b), (c) or (s). Here (Dir-A) has been showed to be well posed by Fabes, Jerison and 
Kenig [H] for small perturbations of (c), using the method of multilinear expansions. 
More recently, the boundary value problems have been studied in the Lp setting 
and for real but non-symmetric matrices in the plane, i.e. n = 1. Here Kenig, 
Koch, Pipher and Toro [20] have obtained solvability of the Dirichlet problem for 
sufficiently large p, and Kenig and Rule ^22j have shown solvability of the Neumann 
and regularity problems for sufficiently small dual exponent p'. 

In the perturbed case A ^ Aq when Aq satisfies (c) or (s), the well posedness of 
(Neu-A), (Reg- A) and (Dir-A) is also proved in [2J by Alfonseca, Auscher, Axels- 
son, Hofmann and Kim. With the further assumption of pointwise resolvent kernel 
bounds, perturbation of case (b) is also implicit in [2]. It is worth comparing the 
present methods to those of [2j. In [2], due to the presence of kernel bounds, the 
solvability of the boundary value problems is meant in the sense of non-tangential 
maximal estimates at the boundary and this follows from the use of layer potentials. 
The first main result in [2] in the unperturbed case (s), is the proof via singular 
integral operator theory of boundedness and invertibility of layer potentials. The 
second main result in [2] is the stability of the simultaneous occurrence of both 
boundedness and invertibility, which hold in the unperturbed cases (c), (s) and (b). 
Solvability then follows. 

Here, we setup a different resolution algorithm (forcing us to introduce some sub- 
stantial material), which consists in solving the first order system (11.21) instead of 
(II. ip . also by a boundary operator method, but acting on the gradient of solutions 
involving a generalised Cauchy operator E^, the goal being to establish boundedness 
of Ea and invertibility of related operators Ea±Na- Boundedness of Ea = sgn(T4) 
is obtained via quadratic estimates of an underlaying first order differential operator 
T4, and the deep fact is here that those quadratic estimates alone are stable under 
perturbations. Stability of invertibility is then easy. The perturbation argument 
requires sophisticated harmonic analysis techniques inspired by the strategy of [2j. 
In particular, the latter uses extensively the technology of the solution of the Kato 
problem for second order operators in [3], whereas we utilise here the work of Ax- 
elsson, Keith and Mcintosh [9], which adapts and extends this technology to first 
order operators of Dirac type. Indeed, we note that our Dirac type operators T4 are 
of the form Ua of ^ in the case of block matrices (b) of Theorem 11.11 But Ta has a 
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more complicated structure when A is not a block matrix and we understand how to 
prove boundedness of at the moment only in the cases specified by Theorem II .![ 
We also note that the present paper, like [Hj, makes no use of kernel bounds and 
only needs L2 off-diagonal bounds for the operators, which always holds. 

The boundary operator method for first order Dirac type operators, used here to 
solve second order boundary value problems, was developed in the thesis of Axels- 
son [6] , which has been published as the four papers [71 El [101 E] ■ It covers operators 
on Lipschitz domains as described above and in Example 11.51 The result in [9] pur- 
sued the program initiated by Auscher, Mcintosh and Nahmod in [1], consisting of 
connecting the Kato problem and the functional calculus of first order differential 
operators of Dirac type. As said, it thus applies to the boundary value problems for 
operators of case (b). What is new here is the setup for full matrices encompassing 
the above. We prove also a sort of meta-theorem (see Theorem 11.31) which roughly 
says that the set of matrices for which the needed quadratic estimates on Ta hold, 
is open. 

We also show that non-tangential maximal estimates hold for our solutions. By 
uniqueness in the class of solutions of (11.11) with non-tangential maximal estimates, 
this implies that our solutions are the same as those in [2] for perturbations of the 
real symmetric and constant cases. The non-tangential maximal estimate here also 
yields an indirect proof of non-tangential limits of solutions of (II. 2p which hold for 
the solutions of (11.11) in [2]. We do not know how to prove this fact directly in 
the framework of this article. Note also that we prove here that the non-tangential 
maximal functions have comparable L2-norms for different values of the parameters 
Co and ci, and that the slightly different non-tangential maximal function used in 
[2] therefore has comparable norm. 

Before turning to the method of proof for Theorem 11.11 we would like to stress 
the importance of the final result that the solution operators g ^ F and u ^ U 
depend Lipschitz continuously on L^o changes of the matrix A around Aq. This is 
an important motivation for considering complex A, as the authors do not know 
any proof of this perturbation result which does not make use of boundedness of 
the operators in a complex neighbourhood of Aq. We also remark that we in fact 
prove that A Sa is holomorphic, from which we deduce Lipschitz continuity as a 
corollary. 

In this paper, we shall mainly focus on the boundary value problems (Neu-A) 
and (Reg-A). The reason is that, assuming the Cauchy operator Ea is bounded, we 
prove in section 12.51 that well posedness follows as 

(Reg-A*) « (Neu^-A) =^ (Dir-A). 

That (Reg-A*) implies (Dir-A) has been proved by Kenig and Pipher fIT\ Theorem 
5.4] in the case of real matrices A. 

1.1. Operators and vector fields. We now explain the basic ideas of the method 
we use for the proof of Theorem ll.il The appropriate Hilbert space on the boundary 
R'^ is 

:= {/ G L2(R"; C"+^) ; curU(/||) = 0}. 

The condition on / means that its tangential part is curl free. Indeed, the trace 
/(x) of a vector field F(t,x) solving (11.21) belongs to due to the second equation 
in (II. 2p The basic picture, building on ideas from 0, is that the Hilbert space splits 
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into two different pairs of complementary subspaces as 

(1.4) ii} = E+n^ © E^n^ = N+n^ © Nxn\ 

We first discuss the splitting into the Hardy type subspaces E^H^, consisting of L2 
boundary traces of vector fields solving f ll.2p in R!^"*"^ respectively. Our main 
work in this paper is to establish boundedness of the projection operators E^ for 
certain A. These projections can be written E^ = ^{I±Ea), where Ea for simple A 
is a singular integral operator of Cauchy type. However, in the general case Ea may 
fail to be a singular integral operator. To handle the projections E^ we make use of 
functional calculus of closed Hilbert space operators, and show that E^ = x±(^a) 
are the spectral projections of an underlaying bisectorial operator Ta in 7i^. The 
functions X±i^) ^^^'g the characteristic functions for the right and left complex half 
planes. To find T4, assume F(t, x) satisfies (11.21) in R""*"^ and solve for the vertical 
derivative 



n 

dtFo = -aoV ( ("oidiFo + diiAF), 

i=l 



dtFi = diFo, i = l,...,n. 

The right hand side defines an operator — T4 in Ti^ which on F{t, x), for fixed t > 0, 
satisfies 

dtF + TaF = 0. 

Concretely, if we identify / = /oCq + /y with (/o, /y)*, where /y is a tangential curl 
free vector field, then 



;i.5) tai 



Aooi((Aoi,, V,) + div.Aiio) Aoo'diVxA,, 

-v.. 





7o" 




A 



where D(T^) = {/ = (/o,/,,)* G ; V./o G L2, div,(A/)|| e L2} and A 
^00 ^0|l 



^110 ^11 



If F{t,x) is a vector field in R""^ satisfying (11.21) . then using this 



operator T4, we can reproduce F provided we know the full trace / = -F|R,n, through 
a Cauchy type reproducing formula F{t,x) = (e~*'^'*'/)(a;). However, in (Neu-A) 
and (Reg- A) only "half" of the trace / is known since the boundary conditions for 
/ are Cq ■ Af = (p and Cq a / = cq a Vip respectively. 

We now turn to the second splitting in (11.41) . which is used to split the boundary 
trace / into the regularity and Neumann data. We define the A-tangential and 
normal subspaces of 'W to be the null spaces of these two operators. 

N+V} := {fen';eo-Af = 0} 

N^n':={fen' ; eoA/ = 0} 

In contrast with the Hardy subspaces, it is straightforward to show that we have a 
topological splitting = N^'W (BN^'W , and therefore that the corresponding pair 
of projections A^^ are bounded. We can now reformulate (Neu-A) and (Reg- A) as 
follows. The Neumann problem (Neu-A) being well posed means that the restricted 
projection 

: E\V} N^V} 
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is an isomorphism, since N^Ti.^ is a complement of the null space of eo-A{-). Similarly 
the regularity problem (Reg-A) being well posed means that the restricted projection 

N+ : E+n^ — > N+n^ 

is an isomorphism, since N^H^ is a complement of the null space of cq a (•). Note 
that what is important here is which subspace projects along, not what subspace 

they project onto. 

We shall also find it convenient to use the operators Ea '■= E^ — E^ = sgn(TA) 
and Na ■— — N^. These are reflection operators, i.e. E\ — I and A^^ = /, and 
we have E^ = \{I ± Ea) and = \{I ± A^^)- 



Example 1.2. Let n = 1 and A = I. Then the space ii} is simply L2(R; C^) and 
the fundamental operator Ta becomes 






_d_- 

dx 




' 




d_ 

dx 













if ~ denotes conjugation with Fourier transform. Furthermore 

sgn(r) 



E '. — Et 





-iH 



iH 




N := Ni 



-1 




where Hf{x) := ^p.v. J ^z^dy. Note that the operator E is contained in the Borcl 
functional calculus of the self-adjoint operator T, and it follows that = 1. On 
the other hand the operator N is outside the Borel functional calculus of T. Indeed, 
the operators b{T) in the Borel functional calculus of T all commute, but we have 
the anticommutation relation 



(1.6) 



EN + NE = 0. 



A consequence of this equation is that the boundary value problems (Ncu-/) and 
(Reg-/) are well posed, or equivalently that A^^ : E^L2 —>■ N~L2 and A^"*^ : E^L2 — > 
N'^L2 are isomorphisms. Consider for example (Reg-/) and assume we want to solve 
A^+/ — g ioT f e E^L2. Applying 4£'+ to the equation gives 

AE+g^{I + E){I + N)f^f + Ef + Nf + ENf^f + f + Nf-Nf^2f, 

so / = 2E~^g and it follows that A^"*" : E^L2 N^L2 is an isomorphism. Having 
solved for / e E'^L2, we can find the solution F{t, x) = C^f{x) in by using the 
Cauchy extension :— e~*l^l£^"^ for t > 0. As a convolution operator, the Fourier 
multiplier has the expression 



{uoeo + uiei) 



and in particular 



_]_ f tuo(y) -{x- y)ui{y) 



t^ + {x- yY 



dy j Co 
{x - y)uo{y) + tui{y) 



F{t,x) = C+f{x) = 2C+{g^e^) = - 

./R 



t2 + (x - y)2 

X - y)eQ + tei)gi{y) 
f + {x- yf 



dy ei. 



dy. 
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For a more general A, even if A is real symmetric, the operator Ta is not self- 
adjoint and proving boundedness of is a highly non trivial problem. Also the 
equation (11.61) fails for general A, and therefore such explicit formulae for the solution 
F{t, x) as in Example 1 1.2 1 are not available. However, to show well posedness of (Neu- 
A) and (Reg-y4) it suffices to show that / ± ^{EaNa + NaEa) = \{Ea± NaY are 
invertible, as explained in To summarise, in order to solve (Neu-A) and (Reg- A) 
we need 

(i) that the Hardy projections E^ are bounded, so that we have a topological 
splitting = EXV} ® E^ii}, and 

(ii) that the restricted projections NX : EXH^ N^ii^ and : E\ii^ 
N^IH} are isomorphisms. 

To prove Theorem ll.il we shall prove the following. 

(i') That Ta satisfies quadratic estimates for all A such that ||A — Ao||oo < ^■ 
From this it will follow that E^ are bounded for all such A and that A i— > E^ 
is continuous (in fact holomorphic) . 
(ii') That (ii) holds for the unperturbed Aq. From (i') and since clearly are 
bounded and depend continuously on A, it then follows from a continuity 
argument that (ii) holds for all A in a neighbourhood of Aq. 

We emphasise that boundedness of the Hardy projections E^ alone does not show 
that (Neu-A) and (Reg-A) are well posed. In our framework, the boundary value 
problems being well posed means that the Hardy space E^n, which exists as a 
closed subspace when the Hardy projections are bounded, is transversal to the A- 
tangential and normal subspaces N^n. A concrete example showing that (Neu-A) 
and (Reg-y4) may fail to be well posed, even though Ea is bounded, is furnished by 
the matrix 

1 fcsgn(a;) 
-k sgn(x) 1 

with parameter /c G R. The corresponding elliptic equation (11.11) in was studied 
by Kenig, Koch, Pipher and Toro in [20j, where they showed that (Dir-A) fails to 
be well posed for certain values of k. Moreover, that (Neu-A) and (Reg-A) also fails 
for some fc, is shown by Kenig and Rule [22] • On the other hand, \\Ea\\ = 1 for all 
G R since according to (ll.Sp 

^ rfc(sgn(x)|:-^sgn(a:)) £ 
-T- 

ax 

is self-adjoint, and therefore E^ are orthogonal projections. 

1.2. Embedding in a Dirac equation. Unfortunately there is a technical problem 
in applying harmonic analysis to the operator Ta in order to prove (i'): the space 'W 
is defined through the non local condition curla;(/||) = 0. This prevents us from using 
multiplication operators, for example when localising with a cutoff f ^ rjf , as these 
does not preserve H^. To avoid this problem we embed C H := L2(R"; AcR""*"^), 
where 

AcR"+^ = A° © © © . . . ® A"+^ 
is the full complex exterior algebra of R"+^, which in particular contains the vectors 
= C""*"^ and the scalars A*^ = C along with all fc- vectors A'^. (We identify k- 
vectors with the dual /c-forms in euclidean space.) In this way we obtain closure, 



A{x) 
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i.e. all operators, including multiplication operators preserve Ti. Furthermore we 
embed the equation fll.2l) in a Dirac type equation 

(1.7) {dt,. + B{xy^dl^B{x))F{t,x) = 0, 

where B, B E Loo(R"; >C(AcR-'^'^^)) are accretive and B\^i = A. Here dt^x and 
d*^ are the exterior and interior derivative operators, as defined in section [2j In 
particular, if F{t, x) : R""*"^ — >■ is a vector field satisfying (11. 2p then 

dl^iBF) = -divt^AF) = 0, 
dt,xF = curlt^a;F = 0. 

Thus (11.71) follows from (11.21) for any choice of the auxiliary function B. In the same 
way as in Section 11.11 we shall solve for the vertical derivative in (ll.7p and obtain an 
operator Tg acting in 7i, such that = Tsly^i- For applications to the Neumann 
and regularity problems it suffices to consider B = / ©A©/©...©/. The stability 
result for the operator Ta we prove in this paper is the following. 

Theorem 1.3. Let Aq be at-independent, complex, accretive coefficient matrix func- 
tion such that Tbq has quadratic estimates in Ti, where Bq = /©Ao ©/©...©/. 
Then there exists eo > depending only on ||Ao||oo; ^^^^^ such that if A E 
Loo(R"; /^(C""*"^)) is t-independent and satisfies \\A — Ao||oo < Co; then Tb has qua- 
dratic estimates in Ti, where B = I®A®I®---®F In particular Ta has quadratic 
estimates in 'H} . 

Thus Ti^ splits into Hardy subspaces, the spectral subspaces of Ta, i.e. each f G 
Ti,^ can be uniquely written f = f~^ + f~, where = F^\-^n and F^{t,x) = 
e^*'"^'*'£'^/(x) satisfies jil.^) in R!^~^^. Moreover, we have equivalence of norms 

ii/ib^iirih + iirih 

and ^ sup,>o \\Fith ^ \¥dtFi,\\\ ^ ||iV,(F±) l^. 

With a more general choice 

5 = fiO © fii © © . . . © 

where B^ G Loo(R"; £(A^)), we also obtain new perturbation results concerning 
boundary value problems and more generally transmission problems for /c-vector 
fields. For fc-vector fields we consider the function space 

n% ■= {/ e L2(R^ c{a')) ; 4(/ii) = = diUB'f)^)} c n, 

where /y and /x denote the tangential and normal parts of /. This is the appropriate 
function space since traces of fc-vector fields F(t,x) in R^."^^ satisfying 

^g^ {dl,{B''ix)F{t,x)) = 0, 

\dt,xF{t,x) = 0, 

belong to H%. Note that for 1-vectors, i.e. vectors, the condition dl{{B^ f) j^) = is 
void and T-i^ = n. For fc- vector fields, we consider the following. 

THE TRANSMISSION PROBLEM (Tr-B'^a^). 

Let B^ = B^{x) G Loo(R"; ^(A^)) be accretive and let a± G C be given jump 
parameters. Given a /c- vector field g G 7i^, find fc- vector fields F^{x) = F~^{t,x) in 
R:;+^ and Fi{x) = F-{t,x) in Rl+^ such that F^^ G Ci(R+; L2(R"; A'^)) and F± 
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satisfies fll.Sp for ±t > 0, and furthermore \imt^±ao = and lim(^o± Ft = 
L2 norm, where the traces satisfy the jump conditions 

Co A («-/+ - = eo^g 

eo^{B\a+r -a-f-)) = eoAB^g). 

The second main resuh of this paper is the following L^o perturbation result of 
the transmission problem (Tr-B^a^). 

Theorem 1.4. Let Bq = Bq{x) G Lao^RP", C{A^)) be t-independent, accretive and 
possibly complex, and assume that Bl is a block matrix, i.e. 

{BQ)st = 0, whenever G (s \ t) U (t \ s) 

and s,t G {0,1, ... ,n} has lengths |s| = |t| = A;. Then there exists e > and C < 00 
depending only on ||i?Q||oo, the accretivity constant n^k and dimension n, such that 
if B^ = B^{x) G Loo(R."", jC{A^)) is t-independent and satisfies 

(1.9) \\B' - 5o'||l^(R") < min (e, C\ia+/ay + 1|), 

then the transmission problem (Tr-B^a^) above is well posed in the sense that for 
every boundary function g G Ti^, there exists unique k -vector fields F^{t,x) with 
properties as in (Tr-B^a^). The solution depends continuously on g with equiv- 
alences of norms \\g\\2 ~ 11/^ + /^Ib ~ ll/^lb + ll/^lb OLnd 

l|/±|h^sup||F±||2^|||t9iF±|||. 

Moreover, the solution operators S^k{g) = F depend Lipschitz continuously on B^ , 
i.e. there exists C < 00 such that 

- ||l2(R")^a' < C\\B2 - -Bf ||L^(Rn) 

when M.9\) holds for B\ and -B|, where \\F\\x denotes one of the two norms above. 

This perturbation theorem for transmission problems for fc-vectors has two im- 
portant corollaries. On one hand it specialises when = 1 to a generalisation of 
Theorem 11.1( b). giving perturbation results for transmission problems across R" 
for the divergence form equation (11.11) . The details of this Neumann-regularity 
transmission problem is stated as (Tr-Aa^) in section l473l 

On the other hand it specialises when either a"*" = or a~ = to a generalisation 
of Theorem 11.1( b). giving perturbation results for boundary value problems for k- 
vectors. Our result for these boundary value problems {Noi-B^) and (Tan-i?'^) is 
given as Corollary 14.171 in Section 14. 3[ 

Example 1.5. In the case B = B, operators of the form dt^x + B^^d^^B appear 
naturally when pulling back the unperturbed Hodge-Dirac operator d + d* with 
a change of variables. As above, consider the region f2 above a Lipschitz graph 
S = {{t,x) ; t = g{x)}. We define the pullback of the field F : Q AR"+i to be 
the field 

{p*F){t,x) ■.= p^{x)F{p{t,x)) 
in R""*"^, where p{t,x) = {t + g{x),x) is the parametrisation of f2, having differential 



12 



PASCAL AUSCHER, ANDREAS AXELSSON, AND STEVE HOFMANN 



acting on vector fields and extended naturally to AR"+^, and p^(x) denotes the 
transposed matrix. From the well known fact that dt^x commutes with p*, we get 
the intertwining relation 

(1.10) K,. + Gdl^G-^WF) = p*{{d,,x + dl,)F), 

where G{x) = {gij{x)) = p^{x)p{x) is the metric for the parametrisation, being real 
symmetric. Solving for the vertical derivative in the equation {dt^x + dtx)-^ — 
in Q, gives us an operator Ds in L2(S; AR""*"^), which is similar to the operator 
(eo - Vg{x))-^{dx + dl) in H. From ffTTOD it follows that 

TG-i(p(xf/(x)) = p(a;f(D2/)(a;). 

It is known that the operator satisfies quadratic estimates, and therefore so does 
Tq-1. For references and further discussion of D^, see [21 Consequence 3.6]. From 
this we get the bounded Clifford-Cauchy singular integral operator 



2 f {g{x)eo + x) - {g{y)eo + y) 



where (T„ is the area of the unit n-sphere in R""*"^ and Eq-i = sgn(Tc-i). For this 
reason, we shall refer to Eb = sgn(TB) as generalised Cauchy integral operators and 
Eb = X±{Tb) as generalised Hardy projection operators, also when B ^ G^^ and 

1.3. Outline and remarks. In section|2]we explain how we use the exterior algebra 
AcR"'*'^ and the exterior and interior derivative operators d and d* . In section [271] we 
introduce the Dirac type operator Tb which extends T4 to the full exterior algebra 
as well as projection operators which extend the A-tangential and normal pro- 
jections from above. Section is concerned with the spectral properties of Tb, 
where we prove that Tb is a bisectorial operator and that the resolvents {XI — Tb)~^ 
has L2 off-diagonal estimates. Section 12.31 surveys the theory of functional calculus 
of bisectorial operators like Tb- In section [231 Lemma [2.491 and 12.55! characterise 
the classes of solutions Ft and Ut respectively to the boundary value problems, and 
are used in particular to prove uniqueness. 

In section [3] we prove (i') quadratic estimates and (ii') invertibility in the unper- 
turbed case B = Bq. The most involved case is when Bq is real symmetric. In order 
to prove the quadratic estimates we use the results from section 12.41 that Tbq leaves 
certain subspaces and Hbq invariant and that therefore it suffices to establish 
quadratic estimates in each subspace separately. 

Section [H contains the core harmonic analysis results of the paper, where we 
establish quadratic estimates for Tb when B k, Bq through a perturbation argument 
based on equation fl4.3p - fl4.8l) . Section [TT] treats the case when Bq is a block matrix 
and makes use of techniques from the solution of the Kato square root problem [3] . 
In Lemma [4.101 we construct a new set of test functions /q which also can be used in 
the Carleson measure estimate in [9] to simplify the proof there. Section 14.21 treats 
matrices of the form B = I®A®I®---®F The key technique is Lemma 14.14! 
where we compare Bq with a corresponding block matrix Bq for which the results 
from section !4.1! applies. The reason why this approximation Bq Bq works is that 
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the normal vector component has additional regularity by Lemma [2. 121 when F 
satisfy the Dirac type equation (11.71) . 

The paper ends with section 14. 3[ where we bring the results together and prove 
Theorems 11.41 11.11 and 11.31 

The reader who is interested only in Theorems 11.11 and 11.31 may ask whether it is 
necessary to work within the full exterior algebra as described in section [T!^ Indeed, 
instead of proving resolvents bounds for through Tb as done in Proposition 12.221 
a direct proof can be based on (II. 5p . To prove a priori estimates < ||/|| when 

(J — irT^ju = /, we multiply this equation with the diagonal matrix ^ 

and obtain the identity 

(1.11) (a 

With an argument as in Lemma [2.21( 11). a priori estimates < ||/|| can be deduced 
from this by proving that the ranges of the two differential systems on the left hand 
side in (11.111) are complementary subspaces in Ti.^, and similarly for the null spaces. 

The problems when working within starts when making use of multiplication 
operators, for example when proving off-diagonal bounds in Proposition 12.251 and in 
the perturbation formula (I4.3p - (l4.8p . but in particular in the proof of Theorem I4.12[ 
It is not clear to the authors how to avoid using at least part of the full exterior 
algebra AcR-'^'^^ at these points. Proving off-diagonal bounds for the resolvents on 
in Proposition 12.251 seems to require having (J + UTb)'^ defined as a bounded 
operator on all vector fields. Also, in (14. 3p - (14. 81) the last factors in each term preserve 
7i^, whereas the middle factor only preserve vector fields, not the curl freeness. 
Therefore we need to have the first factors defined on all vector fields, and moreover 
off-diagonal bounds are needed here. Going back to Proposition 12.251 this seems to 
require having (/ + UTb)'^ defined as a bounded operator also on A^. 

More importantly, to derive quadratic estimates for the perturbed operator T4 
on 'W from (l4.3p -f H78l) . we must assume quadratic estimates for Tb,,, where Bq = 
I(BAq(BI(B---(BI, on all 7i. The reason is that when proving Theorem 14.121 
it does not suffice to consider f E 'U} as the term |||(5t|||op on the right hand side 
is the norm on a bigger space than since 8 does not preserve this space. This 
explains the hypothesis in Theorem 11.31 as well as the necessity of proving quadratic 
estimates on all Ti in section 13.31 

We also point out that for /c-forms, the situation is worse as the spaces Ti^B depend 
on B for k > 2, which makes it impossible to perform perturbation B ^ Bq within 
this space. 
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2. Operator theory and algebra 

In this section we develop the operator theoretic framework we use to prove the 
perturbation theorems stated in the introduction. In particular we introduce our 
basic operator Tb, along with perturbations of the normal and tangential projection 
operators N' and N+. These all act in the Hilbert space H := L2(R"; AcR''+^) 
on the boundary R" = ^R"^^ = 9R"^^. In R"+^ we write the standard ON basis 
as {cq, ei, . . . , e„}, where Cq denote the vertical direction and ei,...,e„ span the 
horizontal hyperplane R". We write the corresponding coordinates as xq, xi, . . . ,Xn 
and we also use the notation t = xq. The corresponding partial derivatives we write 
as di — ^ and dt — do — ^. Our functions f E H take values in the full complex 
exterior algebra over R"+^ 

A - AcR"+^ = © © . . . © A"+^ 

This is a 2""*"^ dimensional linear space with n + 2 pairwise orthogonal subspaces 
A'^ of dimensions ("^^) • With the notation := e^^ a . . . a e^^ if s = {si, . . . , s^} C 
{0,1,..., n} and si < S2 < . . . < s^, the space A*' of homogeneous k-vectors is the 
linear span of {cs ; |s| = k}. In particular, identifying 60 with 1 and the singleton 
set {j} with j, we have A*^ = C and A^ = C""'"^. A general element in A is called a 
multivector and is a direct sum of /c- vectors of different degrees k. 

Definition 2.1. Introduce the sesqui linear scalar product 

\ s t J s 

on A and the bilinear scalar product f ■ g = fs9s- Define the counting function 
a{s,t) := ^{{si,tj) ; Si > tj}, where s = {si}, t = {tj} C {0, 1, . . . ,n}. 

(i) The exterior product f Ag is the complex bilinear product for which 

Cg ACt — (—1)'^^*'*) Csut if s n i = and otherwise zero. 

(ii) The (left) interior product f _i g is the complex bilinear product for which 
(e^ _i Ct, Cu) = {ct, Cg A e,u) for all s,t,u G {0, 1, ... , n}. Explicitly we have 

Cg JCt — (— et\s ii s C t and otherwise zero. 

Example 2.2. The most common situation is when forming a product between a 
vector a — Yl^=o ^j^j ^ ^ /c- vector 

E /.,.....e,,A...Ae.,. 

f)<s\<...<Sk<n 

In this case the interior product a j / is a (A; — 1)- vector, whereas the exterior 
product a a / yields a + l)-vector. Note also that we embed all different spaces A'^ 
of homogeneous /c- vectors as pairwise orthogonal subspaces in the 2"+^ dimensional 
hnear space A. Thus we may add the two products to obtain (the Clifford product) 



a^f + UAfe A^-^ © a'=+^ C a. 
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In the special case when / = 6 is also a vector, i.e. A; = 1, we have 

n 

a = ajbj G A°, 

j=0 

a A 6 = {ttibj — ajbi) Cj a ej G A^, 

0<i<j<n 

where we write = bj. We see that a j 6 coincide with the bilinear scalar product 
a-b. Furthermore, in three dimensions n+1 = 3, the exterior product aAb G can be 
identified with the vector product a x 6 G A^ by using the Hodge star identifications 
e{i,2} ~ eo, -e{o,2} ~ ei and e{o,i} ~ eg. 

The following anticommutativity, associativity and derivation properties of these 
products summarise the fundamental algebra we shall need in this paper. 

Lemma 2.3. If a, 6 G A^ are vectors and f , g, h ^ A, then 

a Ab = —b A a, a aq = 0, 
f A{g Ah) = {f Ag) Ah, f ^ {g ^ h) = {g a f) ^ h 
a j{b A f) = {a ■ b)f — 6 a (a j /) 

We shall also frequently use that if a G A^ is a real vector, then {a_if, g) = (/, aAg). 

Proof. That aAb = —b a a and a a a = is readily seen from Example 12.21 These 
and the associativity f A{g aK) = (/ a (7) a /i are well known properties of the exterior 
product. To see how f j {g jh) = {g a f) jh follows, note first that by linearity it 
suffices to consider the case when /, g and h all are real, i.e. have real coefficients 
in the standard basis {e^}. Under this assumption we pair with an arbitrary w G A 
and use that left interior and exterior multiplication are adjoint operations. We get 

{f ^{g ^h),w) = {g ^h, f aw) = {h,g A{f A w)) 

= (h, {gAf)Aw) = {{g Af)^h,w). 

For the derivation identity, by linearity it suffices to prove 

Ci J {Cj A Cs) + Cj A (Cj J Cs) 

for all s C {0, 1, . . . , n}. This is straightforward to verify from Definition 12. 1[ □ 

Definition 2.4. For a multivector / G A, we write fif := Cq a /, fi*f := Cq j f and 
mf := Co A f + Co ^ f = ifi + fi*)f. We call 

f± = N~f := nn*f = mij*f = fimf, 

f\\ = N~^f ■= l^*f^f = ^l^f = /^*"^/, 
the normal and tangential parts of / respectively. 

Concretely, if / = J2{si,...,s^,} fsu...,sk e^i a . . . a e^^, then its normal part is 

— 5^ fsi, 
{si,...,Sfc}90 
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,--;Sk A . . . A Cg^. 



and its tangential part is 

/II = E 

In particular is normal if G s and tangential if ^ s. Note that both the 
subspace of normal multivectors A^~A and the subspace of tangential multivectors 
iV+A have dimension 2"', although the subspaces N^a'' have different dimensions in 
general. 

Throughout this paper, upper case letters denote A-valued functions F(t,x) = 
Ft{x) in the domain R""^^ or R"'*'^ whereas lower case letters will denote A-valued 
functions / on the boundary R". We use the sesqui linear scalar product (/, (?) := 
J-^„{f{x),g{x)) dx on the Hilbert space TC. 

Definition 2.5. Using the nabla symbols V^; = V = X]j=i ^j'^j ^t,x = 

YTj=o^j'^j R*^"*"^, we define the operators of exterior and interior derivation as 

n n 



n n 



dt,.F = Vt,. A F := ^ e, A djF, dl^F = -Vt,, j F := - ^ e,- j djF. 

j=0 j=0 

We shall also find it convenient to use the operators d := imd and d* = —id*m = 
imd*. 

Remark 2.6. Clearly, in the splitting A = A° © A^ © . . . © A"+^, the operators d 
and d* map d : L2(R"; A^) ^ L2(R"; A^+^) and d* : LiiR"; A^) ^ LiiR"; a''-^) as 
unbounded operators. Moreover, if we further decompose the space of homogeneous 
fc-vectors into its normal and tangential subspaces as 

A = aJ © (A^ © aJ) © (A^ © aJ) © . . . © (A'^ © a;) © Al^\ 

where A^^ := N^a'' and A^ := N^A^, then we see that the operators d and d* map 

d : L2(R"; A^) L2(R"; Af), d : L2(R"; Af) L2(R"; A^+'), 
d* : L2(R"; Af) — ^ L2(R"; A^), d* : L2(R"; A^) — ^ L2(R"; Af-^^). 

Lemma 2.7. We have, on appropriate domains, {dt^xY — i'^txY = d"^ = (d*)^ = 
= (/U*)^ = and the anti commutation relations 

{m, dt^x] = dt = {m, -d*,^}, mn? = {fx, jj*} = I, 

n 

{d, d*} = -A = -J2 dl {d, /x} = {d, /X*} = 0, 

1 

where {A, B} = AB + BA denotes the anticommutator. 

Proof. The proofs are straightforward, using Lemma 12.31 Let us prove that dt = 
{m, —d^^}. For a function F{t,x) we have 

n n n 

Vx,t J (mF) = ^ Ci J {mdiF) = ^ j (cq j diF) + ^ j (cq a diF) . 
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Using the anticommutativity and associative properties, the first sum is 

n n / n ^ 

y^(eo A Ci) J diF = - ^(cj A Co) j diF = -cq j j diF 

\ 

whereas using the derivation property, the second sum is 

n / " \ 

^ Ci J (eo A diF) = doF - cq a j ^ j S^F j . 
\ / 

Adding up we obtain V^,* j (mF) = dtF — m(Vx,t F)- D 
2.1. The basic operators. 

Definition 2.8. Throughout this paper we denote by 5 G Loo(R."; ^{^)) ^ bounded, 
accretive and complex matrix function acting on / G 7^ as f{x) i— >■ B{x)f{x), with 
quantitative bounds ||-B||oo and kb > 0, where is the largest constant such that 

Re{B{x)w,w) > k,b\w\'^, for all w G A, x E R"'. 

We shall also assume that B is of the form i? = 5° © 5^ © i?^ © ... © where 
B'' G Loo(R"; /^(a'^)), so that B preserve the space of fc-vectors. For the matrix 
part B^ acting on vectors, we use the alternative notation A = B^. We also define 
B := mBm. It is not true in general that B preserve fc-vectors. 

Remark 2.9. It would be more optimal to replace the quantitative bounds ||-B||oo 
and Kb > with Kb and fc^, where Kb and ks are the optimal constants such that 

kB\\f\\'<\{BfJ)\<KB\\f\\\ for all fen. 

We consider F{t, x) satisfying the Dirac type equation 

(2.1) {mdt^x + B-^mdl^B)F{t, x) = 
or equivalently 

(2.2) {dt,, + B-^dl,B)F{t,x) = 0. 

In order to solve for the vertical derivative dtF, we note that 
dt,^F = dF + fidtF, dl^F = d*F - ^i*dtF. 
Inserted in (12.11) this yields 

^tF + jv+V-ijv-g M + B-^md*B)F = 0. 
Definition 2.10. Write Mb '■= iV+ — B^^N^ B and define the unbounded operator 
Tb := MB^md + B-^md*B) = -iMj^\d + B-^d*B) = (rf + ^~^rf*E) 

in n with domain D{Tb) := D{d) n B-^D{d*). 

A rather surprising fact is that the most obvious choice for B, namely B = B is 
not the the best, but rather B = mBm. For example, this is the only choice for 
which a Rellich type formula, as in Proposition 13.81 holds on all Ti, in the case of 
Hermitean coefficients B. 

Note that Tg is closely related to operators of the form fl^ = F + B~^r*B, where 
F denotes a first order, homogeneous partial differential operator with constant co- 
efficients such that F^ = 0, which were studied in [9j. Unfortunately, the factor M^^ 
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does not commute with for general B. However, it has other useful commutation 
properties. 

Lemma 2.11. The operator Mb is an isomorphism and 

{B-^N~B)Mb^ = M^^N~, N~M]^^ = Ms\B-^N-B), 
{B-^N+B)M]^^ = M]^^N+, N^M^^ = M^\B-^N+B). 

Proof. Note that Mb = B^^{BN^ — N^B), where the last factor is the diagonal 
matrix {-B^^) © in the splitting H = N-H © N+H, if and denote 
the diagonal blocks of B, and thus BN^ — N~B commutes with N^. Furthermore, 
the diagonal blocks are accretive, so BN^ — N~B and thus Mb is invertible. This 
proves the two equations to the right. To obtain the equations to the left, replace 
B by B^^ and note that 

-A^- + B-^N+B = N+ - B-^N-B. 

□ 

Let us comment on the terminology "Dirac type equation" for (12.11) . Normally 
this denotes a first order differential operator, like dt^x + d^^, whose square acts 
componentwise as the Laplace operator. In our situation, the following holds. 

Lemma 2.12. If F{t,x) satisfies ^2. ij) . then d^ ^Bdt^x{fTiF) = 0. In particular it 
follows that 

diYt,xAVt,xF^'' = 

if F = + (F^'°eo + -F^'") + + . . . + i.e. is the normal component of 

the vector part of F. 

Proof. We use the anticommutation relations mdt^x = dt — dt^x^ and md^^ = —dt — 
dl^m from Lemma [2.71 which shows that (12.11) is equivalent with 

{dt,xm + B~^d;^mB)F{t, x) = 0, 

since dt — B~^dtB = 0. Applying d^,^B to this equation shows that d^^^^Bdt^xi^F) = 
0, and evaluating the scalar A° part shows that diYt,xA'Vt,xF^'^ = since dl^^Bdt^x 
preserves fc-vectors. □ 

The following notions are central in our operator theoretic framework. 

Definition 2.13. Let ?i be a Hilbert space. 

(i) A {topological) splitting of 7i is a decomposition H = Hi (B ^2 into closed 
subspaces Hi and 7^2- In particular, we have ||/i + /2II ~ ||/i|| + II/2II if 
/. e H,. 

(ii) Two bounded operators R'^ and i?~ in H are called complementary projec- 
tions if R+ + R- = I, {R^f = R^ and = 0. 

(iii) A bounded operator i? in is called a reflection operator if R^ = I. 

We note the following connection between these concepts. 

Lemma 2.14. There is a one-to-one correspondence 

R = R+ -R- < — >R^ = 1{I±R)< — >H = R+H © R'H 

between reflection operators in H, complementary projections in H and topological 
splittings ofH. We write R'^H = for the range of the projection i?^. 
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In Definition 12.41 we introduced the complementary projections A^^ associated 
with the sphtting of Ti into the subspaces of tangential and normal multivector 
fields. The corresponding reflection operator is 

:= A^"*" — = n* ji — fifi* = (/i* — fi)m = m{fi — jj*). 

We also introduce 5-perturbed versions of the tangential and normal subspaces 
N-n and N+H as 

B-^N+n := {B-^f ; / G N+H}, 
B-^N-n := {B-^f ; / G N-H}. 

In Definition 12.101 we encountered one of the complementary projections B~^N^B 
and B^^N^B associated with the splitting 

n = B-^N^n © B-^N-n. 

However, more important will be the following complementary projections. 

Definition 2.15. Let and be the complementary projections associated 
with the splitting 

n = B-^N+n®N-n. 

We sometimes use the shorter notation := and iV^ := iV^. Also let A"^ and 
be the complementary projections associated with the splitting 

H = N+n © B-^N-H. 

Let Nb = Nb '■= — and Nb := A"^ — be the associated reflection 
operators. 

With the notation /i^ := B^^^*B, these operators are 

Nb = f^if^ + fJ'*By^ = + 

and similarly for and Nb- We shall prove in section [2?2] that all these operators 
are bounded. It is mainly the operators N^, and Nb that we shall use. 

Definition 2.16. Let (/,c/)b := ((5A^+ - N'B)f,g), for f,g E H. As BN+ - 
N^B = BMb is invertible, (■, ■)b is a duality, i.e. there exists C < oo such that 

K/,^?)b|<C7||/||||^||, 
11/11 <Csup|(/,(7)i.|/||^7||, <Csup|(/, (7)^1/11/11. 

We write T' for the adjoint of an operator T with respect to this duality, i.e. if 
{Tf,g)B = {f,T'g)BioT all f,gen. 

Proposition 2.17. We have adjoint operators (Nb)' = Nb*, (Nb)' = Nb*, N' = N 
and {TbY = -Tb* . 
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Proof. To prove that (Nb)' = Nb*, recall that Nb is the reflection operator for the 
sphtting H = B^^N^Ti © N^Ti and Nb* is the reflection operator for the splitting 
n = N+n © {B*)-^N-n. Thus we need to prove that 

{B'^N^h - N'f2, N+gi + {B*y'N-g2)B 

= {B-'N+f, + N-f,, N+g^ - iB*)-'N-g2)B, 

for all /„ g,, i.e. that {B-^m fi, {B*)-'N- g2)B = = {N' f2, N+gi)B. Use 
Lemma [2.111 to obtain 

{B-'N+fu {B*y'N-g2)B = {{BN+ - N-B)B-'N+fu {B*)-^N-g2) 

{N~MBB-^N+h,g2) = {MBB-'N'N+fug2) = 0. 

A similar calculation shows that {N^ f2, N^gi)B = 0. The proof of (Nb)' = Nb* 
is similar. For the unperturbed reflection operator N, the duality N' = N follows 
directly from the fact the BN^ — N~ B is diagonal in the splitting 7i = N^Ti^N'Ti. 
To prove (Tg)' = —Tb*, we calculate 

(Ts/, g)B = - N-B)M^\md + B~'md*B)f, g) 

= {{Bmd + md*B)f, g) = (/, {B*dm + d*mB*)g) 

= -(/, {N+B* - B*N-){N+B* - B*N-)-\B*md + md*B*)g) 

= -if, (iV+ - {BT'N-BT\rnd + {BT'md* B*)g)B = (/, -TB*g)B, 

where we have used that {B*)-^N+B* - = N+ - {B*)-^N-B*. □ 

Remark 2.18. Our main use of these dualities is for proving surjectivity. Recall 
the following standard technique for proving invertibility of a bounded operator 
T : Til — > '^2- Assume we have at our disposal two pair of dual spaces (7ii,/Ci)i 
and (7i2,A^2)2 and that the adjoint operator is T' : /C2 — >^ JCi. If we can prove a 
priori estimates 

||T/||> 11/11, for all /G Til, 
then T is injective and has closed range. If furthermore T' is injective, in particular 
if it satisfies a priori estimates, then T is surjective and therefore an isomorphism. 

Remark 2.19. We shall also need to restrict dualities to subspaces. Let (7i, /C) be a 
duality, i.e. let (■,■): 7i x /C —> C satisfy the estimates in Definition 12. 161 If Tii C Ti 
is a subspace, then a subspace /Ci C /C is such that (7ii,/Ci) satisfies estimates as 
in Definition 12.161 if and only if /Ci is a complementary subspace to the annihilator 
{g E }C ; {f,g) = 0, for all / G Tii}, in the sense of Definition I2.13( i). 

In particular, if are complementary projections in Ti, the adjoint operators 
(i?^)' are also complementary projections and the duality {H, 7i) restricts to a dual- 
ity {R'^Ti, [R^yTi.), since the annihilator of R'^Ti. is N((i?"'")'), which is complemen- 
tary to {R+yn. 

2.2. Hodge decompositions and resolvent estimates. In this section, we es- 
timate the spectrum of the operator Tb- For this we make use of Hodge type 
decompositions of H as explained below. 

Definition 2.20. By a nilpotent operator F in a Hilbert space Ti, we mean a closed, 
densely defined operator such that R(F) C N(F). In particular F^/ = if / G D(F). 
We say that a nilpotent operator is exact if R(F) = N(F). 
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If r is another nilpotent operator, then we say that T and T are transversal if 
there is a constant c = c(r, F) < 1 such that 

|(/,^7)|<c||/||||(7||, / G R(f) , ^? e R(r) , 

or equivalently if ||/ + g\\ ^ ||/|| + ||5f|| for all / G R(r) and g G R(f). Note 
that any nilpotent operator F is transversal to its adjoint F* with c = 0, since 

R(f) C N(F) = R(F*)^. 

Below we collect the properties of Hodge type splittings which we need in this 
paper. This generalises results from P Proposition 2.2] and (TUl Proposition 3.11]. 

Lemma 2.21. Let F and F be two nilpotent operators in a Hilbert space Ti which 
are exact and transversal with constant c(F,F), and assume also that the adjoints 
F* and T* are transversal. Let B be a bounded, accretive multiplication operator and 
assume that 

max(c(F,f),c(F*,r)) < ^b/\\B\U 

Define operators U := F + F, F^ := B^^TB and TIb '■= T + Tb with domains 
D(n) := D(F) n D(r), D(fB) := B'^Di^) and D{YIb) ■= D{T) n D(Tb) respectively. 

(i) We have a topological splitting 

n = A/(F) © N(Tb), 

so that II /i + /2II ^ II /ill + 11/2 II, /i G A/(F), /2 G N{fB).^ The operator Ub is a 
closed operator with dense domain and range. Furthermore Tb is an exact nilpotent 
operator. 

The complementary Hodge type projections associated with the splitting are 

P^ := Fn^i = U^'Tb = TU-^'Tb, 

p| := f ^n^^ = n^^F = f ^n^'F, 

where we identify a bounded, densely defined operator with its bounded extension to 
Ti. In the special case when B = I we write P^ and P^ for these projections. 

(ii) Let Bi and B2 be two bounded, accretive operators in Ti, with /tBj/||-Bi||oo > 
max(c(F, F), c(F*, F*)), i = 1,2. Then there exists C < 00, depending only on 
||-Bi||oo; i^Bi, c(F,F) and c{T* ,T*) , such that 

||(F + B^'fB2)f\\ > C-'X\\fl for all f G D(F) n B^'Dif), 
^/||n/||>A||/|| /or allfeD{U). 

Proof. Since F^ is conjugated to F, it is an exact nilpotent operator. To prove (i), 
it suffices to prove the estimate 

(2.3) ll/ll + lkll<ll/ + ^?l|, /eN(F), geMiVB). 

Indeed, this shows that N(F) © N(Fb) C H. Furthermore, replacing (F, F, i?) with 
(r*,F*,5^*) shows that N(f*) © N(F^,) C H. Conjugating with B* then shows 
that N((Fb)*) © N(F*) C H. Therefore a duality argument proves the splitting 
H = N(F) © N(fij). From this it follows that 

D{Ub) = {D(Tb) n N(F)) © (D(F) n N(Tb)) 
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and R(nB) = R(r) © R(rB) are dense and that 11^ is closed, as well as the bound- 
edness of the associated projections. 

To prove (EJ]), we use that Bg G N(f) and thus \{f,Bg)\ < c(r, f ) ||/|| and 
estimate 

ll/f < K-^' Re{Bf, /) < K-^' Re + g), /) - {Bg, /)) 

<«:^i||i?|U||/ + ^7|| 11/11 + «:^^c(r,f)||i?^7||||/|| 

<K^^||5|U||/ + ^||||/|| + /c^Mr,f)||5|U(||/ + ^7|| + 11/11)11/11 

</s:^l5|U(l + c(r,f))||/ + (7||||/||+/s:^l5|Uc(r,f)||/||2. 

Solving for ||/||2, this shows that ||/|| < ||/ + g\\ provided ||5||ooc(r, f ) < 1, 
which proves (12.31) . 

To prove (ii) we factorise 

r + 5f ^2 = (p^ + fif ip2)n(p2 + pifia)- 

Here (pi + Sf^P^)-! = P],^ + BiPl^ and (P^ + pi^a)-! = P],^B^' + Pl^. Indeed, 
we have 

(pi+i?r^p')(pk + ^iP|J = p'Pk + p'^iPl. 

+ B^'P^Pl^ + B^^P^BiPl^ = Pip^^ + Efip^SiPl^ 
= (Pi + P2)P^^ + B^\P' + P^)i?iP|, = P],^ + P|, = /, 
{P'+P'B2){P],^B^' + P|J = P'Pl,,^^"^ + P'P|, 

+ P'B2P],^B^' + P'52P|, = P'P|, + P'B^P'b^B^' 

= P'(Pk + PL) + P'^2(Pk + PIJ52-' = P^ + P^ = /. 

A similar calculation shows that P]^^ + i?iP^^ and P^^^B2^ + P^^ also are left 
inverses. Thus P^ + B^^P^ and P^ + P^-B2 are isomorphisms and (ii) follows. □ 

We can now prove the following perturbed normal and tangential splittings of Ti. 

n = N+n © N^n = B'^n+h © N'H 
n = N+n © N^n = n+h © b-^n-h 

To see this, let first T = fi and F = F* = /i* in Lemma l2.21( i). It follows that 
N(r) = N-n, N(F|j) = B-^N+n, P]^ = and P| = iV+. On the other hand, 
choosing F = /i* and F = F* = fi, we see that N(F) = N+H, N(F|j) = B'^N^H, 
P]^ = and P^ = N^. Lemma [2.21( i) thus shows that all the oblique normal 
and tangential projections and from Definition 12. 151 are bounded operators 
on Ti, i.e. we have the stated splittings. 

Next we apply Lemma [2.21( 11) to prove resolvent bounds for the operator Te from 
Definition 12.101 Define closed and open sectors and double sectors in the complex 
plane by 

5"^+ := {z eC ; \ aigz\ < u} U {0}, := U (-5"^+), 

■.= {zeC ■ zj^O, I arg^l < z/}, := U (-S,%). 
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Proposition 2.22. The operator Tb is a closed operator in Ti with dense domain 
and range, for any accretive, complex matrix function B G Loo(R"; £(A)). Further- 
more, Tb is a bisectorial operator with ct{Tb) C S^^, where 

ijj := arccos(KB/(2||-B||oo)) G [7r/3,7r/2). 

and if uj < V < 'njl, then there exists C < oo depending only on u, kb and ||-B||oo 
such that 

\\{X-Tb)-^\\<C/\\1 for all X is,. 

Proof. It is a consequence of Lemma [2.21( i) that the operator Hb = d + B~^d*B is 
a closed operator with dense domain and range. Since Tb = —iM^^IlB with Mb an 
isomorphism, it follows that Tb also is closed with dense domain and range. 

To prove the resolvent estimate, write A = l/iir) where r G S"^^-!^- We first 
prove that < C||/|| if (/ — iTTB)u = /, uniformly for r G S°^^_^. Multiply the 
equation with i{fi — B^^^*B) to obtain 

(r + B-^fB)u = tijj - B-^i2*B)f, 

where T := ifi + rd and F := —ifi* + rd* are nilpotent by Lemma [2. 7[ It suffices to 
prove < \\{r + B~^fB)u\\. 

(i) By orthogonality we have 

||(r + r*)Mf = iirwf + ||r*M||2 = Wi^uf + \T\'^\\duf + 2Re{iiJu,Tdu) 

+ ||/i*-u||^ + Irl^llc?*^!^ + 2'Re{—in*u,Td*u), 

where 

Re(z/iM, rdu) + Re{—ifi*u, rd*u) 

= Re{iTd* fiu, u) + Re(M, iT^d*u) = Re{iT{d*, n}u, u) = 0, 

byLemmaO Thus \\(T+T*)uf = \\{fx+fx*)uf + \Tm{d+d*)uf > \\muf = \\uf. 
In particular F is exact. 

(ii) Next we prove that F and F are transversal, with a bound c < 1 uniformly 
for all r G S^^^-u- -^^ exactness, it suffices to bound (/, g) for f = Tu E R(F) and 

Tg = 0. Furthermore, using the orthogonal Hodge splitting 7i = N(F) © N(F*), we 
may assume that F*m = 0. We get 

(/> g) = (rw, g) = {u, T*g) = {u, {-ifi* + Td*)g) 

= {u,-ifi*g + i{T/T)fj,*g) = 2(r - r)/(2ir)(M,/i*^), 

and thus < 2| sin(argT)|||M||||(y'|| < c||/||||(7|| by (i), where c < kb/||-B||oo since 

7t/2 — V < arcsin(/tB/(2||i?||oo))- A similar argument shows that F* and F* are 
transversal with the same constant c < Kb/||-B||oo uniformly for r G S°^i^__^. 

(iii) To apply Lemma [2.21( ii) it now suffices to prove that ||(F + F)m|| > C^^||?i|| 
uniformly for all r G S°i2_j^. From Lemma [2.21( 1) with B = I we have IKF + F)^!! ^ 

\\Tu\\ + \\Tu\\. Using the Hodge splitting H = N(F) © N(F) it suffices to prove 
11^^11 > \\u\\ for u G N(F), and \\Tv\\ > ||f || for v G N(F). To prove for example the 
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first estimate, write N(r) 9 m = tii + tig G N(r) © N(r*). Then 

> + Ihif -2c||m||||mi|| > (l-c^)||u||^ 

where we have used (i) in the second step and (ii) in the fourth step. 

This proves that < C||/|| if (/ — iTTB)u = f. Since (/ — itTb)' = I — itTb* 
by Proposition 12.171 a duahty argument shows that / — itTb is onto, and the proof 
is complete. □ 

From the uniform boundedness of the resolvents Rf := {I + itTB)~^ for t G R and 
the boundedness of the Hodge projections Pj^ := dH]^^ and := ^^11^^, where 
Ub = d + d*B, d = imd and = B~^d*B, we can now deduce boundedness of 
operators related to Pf = \{R^t + Rf) and Qf = j-iR^t - Rf)- 

Corollary 2.23. The following families of operators are all uniformly hounded for 
t > 0. 





■= {I + itTB)-\ 


tdpf 


= tP^BMBQf, 




Pf 






= zPlMBQf, 




Qf 


■.= tTB{i + en)-\ 


P^M^Hd 


= iQfP'^B^ 




t'npf 


= tTBQf = I-Pf, 


P^M-^H^j, 


= iQfPB^ 




tdQf 


= zP],MB{I-Pf), 


tdP^MgHd 


= P^BMBiPf 


-/)PI, 


t&Qf 


= zPlMB{I-Pf), 


tdPfM^Hd*B 


= P^MB(Pf 


-/)Pb, 


fM^Hd 


= ^{I - )P|, 


tdsPfMsHd 


= PlMBiPf 


- I)PI, 


M^Hd^ 


= z(/ - P,^)P^, 


td^P.^M^Hd*^ 


= PlMBiPf 





These families of operators are not only bounded, but have L2 off-diagonal bounds 
in the following sense. 

Definition 2.24. Let {Ut)t>o be a family of operators on Ti, and let M > 0. We say 

that {Ut)t>o has L2 off-diagonal bounds (with exponent M) if there exists Cm < 00 
such that 

\\UJ\\L,^E)<CM{dist{E,F)/t)-''\\f\\ 

whenever E,F G HJ^ and supp / C F. Here (x) := 1 + and dist (i?, P) := 
inf{|a; — y\ : x & E,y & F}. We write ||t/t||off,A/ for the smallest constant Cm- 
The exact value of M is normally not important and we write ||t/t||ofr, where it is 
understood that M is chosen sufficiently large but fixed. 

Proposition 2.25. All the operator families from Corollary \2. 23\ has Lo off-diagonal 
hounds for all exponents M > 0. 

Proof. First consider the resolvents Rf = (/ + HTb)^^- As we already have proved 
uniform bounds for Rf , it suffices to prove 

Ul + ttTB)-'f\\L,(E) < CM(|t|/dist (P,P))*'||/|| 

for \t\ < dist {E,F). We prove this by induction on M as in [9l Proposition 5.2]. 
Let ?7 : R" [0, 1] be a bump function such that tiIe = 1, supp 77 <Z E := {x E 
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R" ; dist(x,E) < dist(x,F)} and \\Vri\\^ < l/dist(E,F) ^ l/dist(E,F). Since 
the commutator is 

[rjI,R^]=tRfM-\[r^I,d]+B-'[r^I,cf]B)R^, 
where ||[?7/,rf]||, < IIVt^Hoo, we get 

WRffh^iE) < WvRffW = \\[vI,Rm < |t|l|Vr7||oo||i?f/IL,(^), 

where we used that rif = 0. By induction, this proves the off-diagonal bounds for 
R^. From this, off-diagonal bounds for Pf, Qf and I — Pf also follows immediately. 
Next we consider tdP^^ and use Lemma F2.21( i) to obtain 

(2.4) \\tdPffh,,^E) < WvtdPffW < \\[vI,td]Pff\\ + WtdvPffW 

<\t\\\Vv\\oo\\Pffh,iE) + \\tTBVPFf\\ 
< |t|||Vr/|U||Pf /ll^^(^) + \\[7^,tTB]Pff\\ + WTsPffW 

<|i|||Vr/||oo||Pf/IL,(^) + ||gf/IL,(^). 

This and the corresponding calculation with d replaced by d*^ proves the off-diagonal 
bounds for tdPf and td^Pt^. From this the result for PfM^^Hd and P^^M^Hd^ 
follows immediately with a duality argument. Indeed, (M^^d)' = M^ld*^, and 
(M^^d*^)' = M^]d is proved similarly to = —Tb* in Proposition 12.171 

The proof for tdQf , td*^Q^, QfM^^td and QfM^^td*^ is similar, replacing Pt and 
Qt with Qt and I ~ Pt- Finally, the last four estimates follows from a computation 
like (12.41) . for example replacing Pt and Qt with PfM^^td and Q^M^^td proves the 
estimate for tdPfM^^^td. □ 

We finish this section with a lemma to be used in Section |H This lemma is proved 
with an argument similar to that in [T51 Lemma 2.3]. For completeness, we include 
a short proof. 

Lemma 2.26. Assume that {Ut)t>o (^i^'d {Vt)t>o both have L2 off-diagonal bounds 
with exponent M. Then {UtVt)t>o has L2 off-diagonal bounds with exponent M and 

\\UtVt\\off,M < 2^^^"^||t/(||oj(f,J\/||Vt||ojff,Af- 

Proof. By Definition 12.241 we need to prove that 

l|f/*^*/llL.(i.) < (dist (E,F)A)-^||/1| 

whenever supp / C F. To this end let G := {x G R*^ ; dist {x,E) < p/2}, where 
p := dist {E,F). We get 

WUtVtfh^iE) < \\UtiXGVtf)\\L,iE) + \\UtiXW^\GVtf)\\L,iE) 

< c^'WVtfKiG) + cS(p/2t)-"1i^*/l|L.(RnG) 

< C^^C^Hp/2t)-'^'\\f\\+C'^Hp/^r''C^'\\f\\ < 2^^+^C^'C]l(p/t)-^||/||. 

□ 
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2.3. Quadratic estimates: generalities. In Proposition 12.221 we proved the spec- 
tral estimate criTB) C for some angle ui < 7r/2 with bounds on the resolvent out- 
side S^. In this section we survey some general facts about the functional calculus 
of the operator Tb- For a further background and discussion of these matters we 
refer to P[T]. 

Definition 2.27. For uo < v < 7r/2, we define the following classes of holomorphic 
functions / G H{S°) on the open double sector S°. 

^{S:) ■.= {i,eH{S:) ■ 1^(^)1 <Cmin(|^|M^|-^),^G^°, 

for some s > 0, C < oo}, 
Hoo{Sl) := {b e H{S°) ; \b{z)\ <C,ze S°, for some C < oo}, 

FiS:) := {w e His:) ; |^(^)| < Cmaxd^^, kl'^,^ e S^, 

for some s < oo, C < oo}. 

Thus ^(^°) C H^iS^) c C H{St). 

For ip G ^E'(5'°), we define a bounded operator ipiTs) through the Dunford func- 
tional calculus 

(2.5) ^(Tb) ■■=7^ [ ^(A)(A/ - Tsr'dX 

27rz 

where 7 is the unbounded contour {±re^*^ ; r > 0}, u < 9 < u, parametrised 
counterclockwise around 5*0;. The decay estimate on ip and the resolvent bounds of 
Proposition 12.221 guarantee that < 00. 

For general w G F{S:) we define 

where k is an integer larger than s if 1^(^)1 < C maxd^l^*, and q{z) := z{l + 

z^)~^ and := qiTB). This yields a closed, densely defined operator w(Tb) in H. 
Furthermore, we have 



XiWi{Tb) + X2W2{Tb) = (XlWi + \2W2){Tb), 
(2.6) Wi{Tb)w2{Tb) = {WiW2)iTB), 

for all and W2 G ^(5*°). Here T = S means that the graph G(T) is dense in the 
graph G{S). 

The functional calculus w 1— > w{Tb) have the following convergence properties as 
proved in [Tj. 

Lemma 2.28. If bk G Hoo{Sl) is a sequence uniformly hounded on 5° which con- 
verges to b uniformly on compact subsets, and if bk{TB) are uniformly hounded op- 
erators, then 

bkiTB)f ^b{TB)f, for all fen, 
and \\b{TB)\\ < limsup^ \\bk{TB)\\. 

Definition 2.29. The following operators in the functional calculus are of special 
importance to us. 

(1) q^(z) = q{tz) := tz{l + t^^^^-i ^ ^(5°), which give the operator Qf. 

(2) l^;!"* := {z^Y^'^ G F{Sl), which give the operator \TbY. Note that \z\ does not 
denote absolute value here, but 2 1— > |z| is holomorphic on S^. 
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(3) e G Hoq{S"), which give the operator e 

(4) The characteristic functions 

±, . _ f 1 ii±Rez > 
^ ~ \ if ±Re2 < 

which give the generahsed Hardy projections := x^(Tb). 

(5) The signum function 

sgn(z) = x'^iz) - X~iz) 
which give the generahsed Cauchy integral Eb '■= sgn(rB). 

The main work in this paper is to prove the boundedness the projections E^. As 
in Lemma 12.141 if these are bounded then they correspond to a sphtting 

H = E+n © E^H 

of H into the Hardy subspaces E^H associated with the equation (12. ip . That the 
projections are bounded is also equivalent with having a bounded reflection operator 
Eb. 

Definition 2.30. For a function F{t,x) defined in R!^"*"^ we write 

and for short =: |||-^|||- When F{t,x) = (Gj/)(a;) for some family of operators 

(0t)i>O5 we use the notation 

|||0i|||op := sup |||6i/|||. 
11/11=1 

Our main goal in this paper will be to prove quadratic estimates of the form 

(2.7) fwQfffj^wfr, 

for certain coefficients B. We recall the following two basic results concerning qua- 
dratic estimates which are proved by Schur estimates. For details we refer to [T]. 

Proposition 2.31. Let ip G '^{SD he non vanishing on both S^^ and S°_, and 
define ipt{z) := ip(tz). Then there exists < C < oo such that 

C-'\\\Qff\\\ < \mTB)f\\\ < C\\\Qff\\\. 

Proposition 2.32. If Tb satisfies quadratic estimates 1^2. then Tb has hounded 
Hoo{S°) functional calculus, i.e. 

\\b{TB)\\<\\b\U forallbeH^iS:). 
Thus Hoo{S°) 3 b I — > b(TB) G C{T-C) is a continuous homomorphism. 

Before proving quadratic estimates for Tb for certain B in Sections [3] and IH we 
introduce a dense subspace on which the operator 6(Tb) is defined for any b G 
Hoo{S°). 

Definition 2.33. Let Vb be the dense linear subspace 

VB:=[jm\TBnnR{\TBn)cn. 

s>0 
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We see that D(|Tb|'^) fl R(|Tb|'^) increases when s decreases. The density of 
D(|Te|) n R(|Tb|), and therefore of Vb, follows from the fact that 

as (a, (3) — i> (0, oo), for all f ^H. 

Moreover, if 6 G Hoo{S°) and f E Vb then b(TB)f eVbCH. To see this, write 

b{TB)f = m{TBmTB)-'f), 

where tpiz)-^ := (1 + \z\')/\z\'/^ if / G D{\Tb\') n R{\Tb\'). Then ip{TB)-^f G 
H and {b'ip)(TB) is bounded since hip G "^^{3°). Furthermore, if s' < s/2 then 
|Tb|^'(V'(Tb)-7) e n and |Tb|-^'(^(Tb)-7) e 7^, so 6(Tb)/ G D{\TB\'')nR{\TB\''). 

Lemma 2.34. We have an algebraic splitting 

Vb = E^Vb + E^Vb, 

E^Vb(~^E]^Vb = {0}, and f G E^Vb is in one-to-one correspondence with F(t,x) = 
Ft{x) := (e^*l^sl/)(x) m Rl+\ and 

hm Ft = f, lim Ft = 0, 

G D(Tb) = D(rf) n D{d*B), dtFt = -TBFt G L2(R"), ±t > 0. 

Proof. That each f eVb can be uniquely written f = + f~ , where G E^Vb, 
follows from (12.61) . To verify the properties of Ft, it suffices to consider the case 
/ G E^Vb as the case / G E^Vb is similar. Since ze"*'^' G ^(5'°) it follows that 
Ft G D(T8). Moreover, since i(e~(*+'*)l^l — e~*l^l) — > — |z|e~*l^l uniformly on S^, 
it follows from Lemma 12.281 that dtFt = — |Tb|Fj and since Ft G E^Vb we have 
|7s|-Pt = TBFt. 

To prove the limits, assume that / G D(|Tb|*) fl R(|Tb|*) for some < s < 1. 
Writing / = \Tb\~'^u, we see that 

Ft- f = t'ilj{tTB)u, where tpiz) = (e-l"' - l)/\z\\ 

Similarly with / = \Tb\''^v, we get 

Ft = t-'ip{tTB)v, where ■^{z) = |z|'e"l^l. 

Since in both cases tp{tTB) are uniformly bounded in t, using a direct norm estimate 
in (12.51) . it follows that both limits are as t — >^ and t — > cx) respectively. In 
particular, / = limj^o Et is uniquely determined by F. □ 

We now further discuss the quadratic estimates (12. 7p . First note the following 
consequence of the duality = — Tg* from Proposition 12.171 Again we refer to pLj 
for further details. 

Lemma 2.35. // |||gf /||| < ||/|| for all f eU, then \\\Qf f\\\ > \\f\\ for all f eU. 

In Section 13.31 we shall use the following Hardy space reduction of the quadratic 
estimate. This is a technique due to Coifman, Jones and Semmes fl], and adapted 
to the setting of functional calculus by Mcintosh and Qian [23| Theorem 5.2]. 
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Proposition 2.36. Assume that we have reverse quadratic estimates in E^,Vb* and 

\\g\\<\\\tdtGt\U, geE^VB*, 
where Gt = e^^^^'^'^g. Then 

IIIQf/lll< 11/11, f^n. 

In particular, if we have reverse quadratic estimates in both Hardy spaces for both 
operators Tb and Tb- , then \\\Qff\\\ ^ ||/|| ^ / e 7^, and thus Tb and 

Tb" have bounded Hoo{S°) functional calculus for all uj < u < n/2. 

Proof. Assume first that 

\\9\\<\¥dtGt\\\, geE+,VB*, 

and write ipt{z) '■= tze^*^^^ so that tdtGt = —iptiTB*)g by Lemma f2.'d4\ Let f E H 
and define qf{z) := X^Qtiz) so that qf(TB) = E]^Qf . To prove that 

\htiTB)f\\\<\\fl fen, 

it suffices to bound ||% y uniformly for a > and (3 < oo. To this end, 

define the auxiliary functions 

ht := hi{TB)fr f ) ^ {N^B* - B*N-)-\;{TB)f, 

9--=- [%tiTB')htf, 

J a 

so that < C and g G E^,Vb*, and calculate 



13 \ 1/2 

IIC(T^B)/f f j = {qr{TB)f,k)Bf = {f,g) 

9\\ < 11/11 \\mTB^)9\\\ 



B 



< II fll llnll < 



(f(/;„(.Mn,...„„,.,,f)^f) 



1/2 



,dt^ ^^"^ 



< 



< 11/11 ^ IM^jj < 

In the second equality we have used that q^iTB)' = qt{—TB*) = —q^(TB*), in the 
second estimate we used the hypothesis and the second last estimate is a Schur 
estimate. We here use that \\{4's<lt')(TB*)\\ ^ vi't/s), where ri{x) := min(x'^, x"'') for 
some s > 0. 

With a similar argument |||g/'(Te)/||| < ||/||, / G H, follows from the reverse 
quadratic estimate for g G -E^*Vb*- If both reverse estimates holds for B*, then 

IIIQf /III < IIIC(Ti.)/lll + 111%-^ra/lll < 11/11, / G H, 

and if the same holds for B and B* interchanged, then Lemma 12.351 proves that 
IIIQf/lll ~ 11/11, / ^ aiid by Proposition 12.321 this proves that Tb has bounded 
Hoo{S°) functional calculus, and similarly for Tb*. □ 



30 



PASCAL AUSCHER, ANDREAS AXELSSON, AND STEVE HOFMANN 



We end this section with a discussion of the holomorphic perturbation theory for 
the functional calculus of Tg. 

Definition 2.37. Let 2; ^— ^ f/^ G C{X, y) be an operator valued function defined on 
an open subset D C C of the complex plane. We say that Uz is holomorphic if for 
all z E D there exists an operator U'^ e C{X,y) such that 

Lemma 2.38. Let z ^ Uz C{X, y) be an operator valued function defined on an 
open subset D C C of the complex plane. Then the following are equivalent. 

(i) z Uz is holomorphic. 

(ii) The scalar function h{z) = {Uzf,g) is a holomorphic function for all f E X 
and g G y* , where X G X and y* C 3^* are dense, and \\Uz\\ is locally 
bounded. 

In particular, if z ^ are holomorphic on D for k = 1,2, . . . and 

{U^f,g) = hk{z) — . h{z) = {Uzf,g), for all z e D, f e X , g e T, 

and \\Uz\\ < 00 for each compact subset K C D, then z ^ Uz is holo- 

morphic on D. 

Proof. For the equivalence between (i) and (ii), see Kato Theorem III 3.12]. 
To prove the convergence result, it suffices to show that h{z) is holomorphic on D. 
That this is true follows from an application of the dominated convergence theorem 
in the Cauchy integral formula for hk{z). □ 

Below, we shall assume that 2; !—> 5^ is a given holomorphic matrix valued function 
defined on an open subset D (Z C such that Bz is a multiplication operator as in 
Definition 12.81 for each z E D, and that ujn := sup^g^, arccos(KB^/(2||i?^||oo)) < 7r/2. 
Let uj£) <v< 7r/2. 

Lemma 2.39. Forr G S^i2_y, the operator valued function D 3 z 1 — > {I — itTb^)^^ 
is holomorphic. 

Proof. Similar to the proof of Proposition 12.221 we have 

(/ - itTbX' = {BzT + fBz)-H{Bz^l - fi*Bz), 

where T = ifi + rd and F = —ifi + rd*. It is clear that the multiplication operator 
Bzfi — ^J*Bz depends holomorphically on z, so it suffices to show that z ^ {BzT + 
TBz)^^ is holomorphic. Let z E D, write B := Bz and B^j := Bz+w and calculate 

(2.8) {B^V + VB.^)-^ -{BV + VB)-^ 

= -(5^F + VB^)-\B^ - B) (F(5F + VB)-^) 

- {{Bj: + TB^Y^f) {B^ - B){BT + fB)-\ 

The first and last factors in both terms on the right hand side are uniformly bounded 
by Lemma [2.21( i) so we deduce continuity of w ^-^ [B^T + TBw)"^. Furthermore, 
multiplying equation (12. 8p from the right with F, the first term on the right vanishes 
as 

F(SF + fB)-^f = {BY + fB)-^f^ = 0, 
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and we deduce continuity of w i— » {B^^T + TB^) ^T. Thus, dividing equation fl2.8p 
by w and letting w — >• we see that the hmit exists and equals 

-{BY + fB)-^B'{T{BT + TB)-^) - {{BY + fBy^f)B'{BT + fB)-\ 

□ 

Lemma 2.40. Ifi'E "^{S"), then D 3 z \ — > 'ipiTsJ is holomorphic. 

Proof. Let 7 be the unbounded contour {±re^*^ ; r > 0}, ud < 9 < u, parametrised 
counterclockwise around 5"^^^ . By inspection of the proof of Lemma 12.391 we have 

sup |A| ||i((A - Ts^^J-' - (A - TsX') - d.{X - TsX'W ^ 0' 
as w — >• 0. Thus 

Since J^mX)\ If I < 00. □ 

Lemma 2.41. Assume thatTs^ satisfy quadratic estimates |||Qf^/||| ~ ||/||, f eH, 
locally uniformly for z E D. If h E Hoo{S°), ip E '^{S°) and Q < a < {3 < 00, then 
the following operators depend holomorphically on z E D. 

(i) u{x) ^ ib{TB^)u){x) -.n^n 

(ii) u{x) I — y v{t,x) = {b{tTB,)u){x) : H — > L2(R" x A) 

(iii) u{x) ^ v{t,x) = {MTb»{x) : n L2{Kl+\ A) 

Proof, (i) Take a uniformly bounded sequence ipkiz) E "^{S") which converges to 
b{z) E Hoo{S") uniformly on compact subsets of S°. Lemma [2.381 then applies with 
= ^kiTB^) and = b{TB,), using Lemma [23Q] and Lemma [223 

(ii) It suffices by Lemma 12.381 to show that h{z) = ht{z)dt is holomorphic, 
where ht{z) = {b{tTBjf,Gt), for all f{x) E H and G{t,x) E C^(W x (a,/5);A). 
That ht{z) is holomorphic for each t is clear from (i), and for h{z) this follows from 
an application of the Fubini theorem to the Cauchy integral formula for ht{z). 

(iii) Consider the truncations : u{x) 1-^ v{t,x) = Xk(t){ipt(TB^)u){x), where Xk 
denotes the characteristic function of the interval {l/k,k). It is clear from (ii) that 

is holomorphic, and letting — 00 we deduce from Lemma [2.381 that iptiTB^) is 
holomorphic. □ 

Proposition 2.42. Assume that Tb satisfies quadratic estimates |||<5f/||| ~ ||/||, 
f E Ti, locally uniformly for all B such that \\B — -BqHoo < cind let ip E "^^{3°). 
Then we have the Lipschitz estimates 

||6(TbJ - b{TB,)\\ < C\\b\U\B2 - 5i|U, b E H^{S:), 
\\mTB,)-MTBj\\\<C\\B2-B^\U 

when \\Bi - -Boll < e/2, i = 1,2. 

Proof. Let B{z) := Bi + z{B2 — i?i)/||i?2 — -Bi||oo, so that B{z) is holomorphic in 
a neighbourhood of the interval [0, \\B2 — -Bi||oo]- In this neighbourhood, we have 
bounds ||&(Tb(2))|| ^ ll^lloo by Proposition 12.321 and holomorphic dependence on z by 
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Lemma [!j.41( i). Schwarz' lemma now applies and proves that || ^6(Tb(2)) || < ||&||oo 
for all z G [0, II-B2 — -Bi||oo]- This shows that 

||KTbJ-6(TbJ||< '''"|||6(Ti.w)||rft<C7||6|U||52-i?i||oo. 

^0 

The proof of the Lipschitz estimate for ^/'j (Tg) similarly follows from Lemma r2.41( iii). 

□ 

2.4. Decoupling of the Dirac equation. In section [?!T] we introduced the Dirac 
type equation 

(2.9) {mdt,^ + B-^mdl^B)F = 

satisfied by functions F{t, x) : R^"*"^ A. Of particular interest is when both terms 
vanish, i.e. when 

idt^F = , , , idF = -fidtF 

(2.10) \ ' , . or equivalently when \ , , „ . 
^ ^ \dl^{BF) = ^ \d*{BF) = fi*BdtF 

Consider a solution F{t,x) = e^*!^^!/ to (ES]) in R;^+^ as in Lemma 12.341 where 
/,Tb/ G E^Vb- Using dtF = -TbF and Lemma EZH we get 

{mdt,^F)\t=o = m{d - iiTb)! = {md - N+M]^\md + B-^md*B))f 
= M-\{Mb - B-^N+B)md - B-^N+md*B)f 
= M^\-N-md - B-^N+md*B)f = M^\dfi + B^^d*fi*B)f. 

Thus we see that dt,^F = = dl^{BF) at t = if and only if dfif = = d*^*Bf. 
We may also rewrite equation (12. 9p as 

(2.11) {dt,xm + B-^dl^mB)F{t, x) = 0, 

by using mdt^x = dt — dt,x^ and md^^ = —dt — d^^m from Lemma \T7\ Consider the 
case when both terms vanish, i.e. when 

(2.12) J ^ , ^ or equivalently when l^ ^^^^ 

^ ' \dl^{niBF) = \d*{BF) = -^iBdtF 

For F{t,x) = e^*'"^^'/ solving (12.91) we also have 

dt,x{mF)\t=o = dtF\t=o - mdt,xF\t=o = M^\dm + B-^d*mB)f 

- M]^\dii + B-^d*fi*B)f = M^^dfi* + B-^d*fiB)f. 

Thus we see that .^(mF) = = d^^imBF) at t = if and only if d^* f = = 
d*^iBf. 

Definition 2.43. Introduce the closed, densely defined operators 

Tb := -Ms\di2* + B-^d*i^B), Tb := -M]^\d^ + B~^d*^i*B), 
with domains 

□ (fs) := {/ G 7^ ; /i7 G D(d), ^5/ G D(rf*)}, 

□ (Tb) ■.= {fEn;fxfE Did), fx*Bf G D(rf*)}, 
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respectively, and define the closed subspaces 

'HB:={fen; difif) = = d*{fi*Bf)}, 

nB:={fen ; d{fi*f) = = d*{fiBf)}. 
Proposition 2.44. We have a topological splitting ofTi into closed subspaces 

Ti = 1H.B © 'Hb- 

Furthermore we have Tb = Tb + Tb with D{Tb) = D(Tb) H D{Tb), and 

Rifs) = nB = a/(Tb), R{fB) = nB = a/(Tb). 

Thus, if we identify Tb with its restriction to Hb and Tb with its restriction to 
Hb, then these are bisectorial operators with spectral and resolvent estimates as in 
Proposition 

Proof Clearly D{fB) n D{fB) = D{d) f] B-^D{d*) = D{Tb) and fB + fB = Tb- 
Furthermore Lemma [2.21( i) shows that N(Te) = Hb and N(Tb) =Hb- 

(1) To show R(Tb) C Hb, let / G D(Tb) and use Lemma [2.111 to get 

d^i{fBf) = -d^iMg^{d^i* + B-^d*i2B)f = -dmN^Mg^{d^i* + B-^d*^B)f 
= -dmMg^B-^N+B{dfi* + B-^d*^iB)f = -dmMj^\B~^N+ B - N-)dfi* f 

= —dmdn*f = md'^fj* f = 0, 

and similarly 

d*^i*B{TBf) = -d*mB{B-^N-B)M-\d^* + B-^d*ijB)f 

= -d*mBMj^'^N^{d^i* + B-^d*fiB)f = -d*mBM]^\N- - B-^N+B)B-^d*fiBf 

= d*md*^iBf = -m{d*fnBf = 0. 

A similar calculation shows that R(Tb) gHb- 

(2) Next we show that for all / G D(Tb), we have 

\\TBf\\ ^ WTBfW + WTBfl 

Using that Mb and B are isomorphisms, Lemma [2.21( 1) with T = d and T = d* and 
orthogonality + = \\g\\^, we obtain 

WTBfW + WTBfW ^ {\\f^*df\\ + \\fid*Bf\\) + (ll/id/ll + \\fi*d*Bf\\) 

^\\df\\ + \\d*Bf\\^\\TBf\\. 

(3) Clearly HBnHB = {0} and R(Tb) C R(Tb) + RITb). Taking closures, using 
that Tb has dense range in H and using (2) yields 

H = R(fB) © R(tB). 

Thus from (1) it follows that R(Tb) = Hb and R(Tb) = Hb and that H = Hb ® 
Hb. □ 

It follows that Tb is diagonal in the splitting H = Hb®Hb- We shall now further 
decompose the subspace Hb into the subspaces of homogeneous /c- vector field A*^, 
which also are preserved by Tg. The same decomposition can be made for Hb, but 
this is not useful since Tb does not preserve these subspaces. 
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Definition 2.45. Let := i:2(R"; C{a'')) and := H'' H Kb- 

Lemma 2.46. The operator Tb preserve all suhspaces 'H'b, Ti-B, 'H^, . . ., IhC^^ and 
"Hb (but not the suhspaces Ti^ ) and we have a splitting 

nB = nl®n],®nl® fii^\ 

Furthermore, for the operators Nb, Nb and N , we have mapping properties 
Nb-.Hb — ' ^1, Nb-.Hb — > Hb and N -.V} — > 

Proof. As d/i*, B and all preserve it follows that Tb(7t;|) = Tb(7t;^) C 
7i^. To show the splitting of Hb it suffices to prove 

Kb cn% (bHb ®nl ® . . . ® nB^\ 

To this end, let / G Hb and write fk for the A'' part of /. Since d/if = = d*fi*Bf 
we get 

= {diJ,f)k+2 = d{{fif)k+i) = dfiifk), 

= {d*fi*Bf)k^2 = d*{{fi*Bf)k^i) = d*fi*{{Bf)k) = d*fi*B{fk), 

for all k. Thus fk G H'^^ and / G 0^ Ti^. 

To prove the mapping properties for Nb, note that if / = /i + /2 in the splitting 
n = B-^N+H © N-n, then f e Kb H and only if G D{d) and fiBfi G D(d*), 
according to Definition [2331 Clearly /i - /2 = NbU) e Hb H f e Hb- Since Nb 
preserves Ti.^, the desired mapping property follows. The proofs for Nb and are 
similar. □ 

Lemma 2.47. With {■,-)b denoting the duality from Definition \2.1b\ we have dual 
operators 

iTBy = ~fB*, {TB)'=~fB*, 

and restricted dualities {'Hb,'Hb*) b, {'Hb,'Hb*) b and {'H^,'H^*)b for all k. In the 
case k = 1, we shall write the duality as ((H},'H})a- 

Proof. The proofs of = —Tb* and = —Tb* are similar to that of = 
—Tb* in Proposition 12.171 From this we get that the annihilator oi I-Lb = R(7b) 
is N (Tb* ) = 'Hb* which is a complement of TiB*- Thus we see from Remark 12.191 
that {'Hb,'Hb*)b is a duality. The proof of the duality {Hb^H-b*) b is similar. Also, 
since BN^ — N~ B in the definition of {■,-)b preserves Ti.'' we also have dualities 
(7t;|,7t!|.)b for all fc. □ 

Remark 2.48. • The subspace of vector fields with a curl free tangential part 

nl = {f eL2{K';A') ; rf(eo a /) = 0} 

is independent of B and coincides with the space from section 11.11 Fur- 
thermore, the operator Ta there coincides with Tb |^i = Tb 
• The dense subspace Vb C H splits = © Vb with Proposition 12.441 
where Vb '■= Vb HTYb C Hb and Vb '■= Vb HTtIb C Hb are dense suhspaces, 
and we can further decompose Vb into homogeneous /c-vector fields 

= Vb © Vb © . . . © Vb+\ 

where V| := Vb n7i:|. 
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• Furthermore, all these dense subspaces V^, Vb, Vs of 7i%, Hb, 'Hb splits 
algebraically into Hardy spaces similar to Lemma 12.341 e.g. = E^V% + 
E^V%. Here / G E%Vb if and only if F{t,x) satisfies fl2inD . and / G E%Vb 
if and only if F{t,x) satisfies (12.121) . 

2.5. Operator equations and estimates for solutions. Our objective in this 
section is to set up our boundary operator method for solving the boundary value 
problems (Neu-A), (Reg- A), (Neu~^-A) and (Dir-A) as well as the transmission prob- 
lem (Tr-B'^a^). Under the assumption that Tg has quadratic estimates, which is 
made throughout this section, we first show that solutions F{t, x) are determined by 
their traces / through the reproducing Cauchy type formula Ft = e~*'"^^'/. Recall 
from Lemma [2.461 that both operators Eb and Nb preserve all subspaces 7i^. We 
shall write EBk = -Esl^fe = sgn(TB|^|) and A^'^fe = A^bI-^^i for the restrictions. In 
particular, we write A = when k = 1. 

Lemma 2.49. Assume that Tb satisfies quadratic estimates, let f G Ti.'^ and let 
(0, oo) 3 t \—>- Ft{x) = F{t, x) G Ti.^ be a family of functions. Then the following are 
equivalent. 

(i) f eEln^ andFt = e-'\''^\f. 

(ii) Ft G C^(R+;7i'^) and satisfies the equations 

(dl^{B{x)Fit,x)) = 0, 
\dt,.F{t,x) = 0, 

and have L2 limits lim(_,o+ Et = f and lim(_>oo Ft = 0. 
In fact, such Ft belong to C-' (R+;7Y^) for all j > 1 and are in one-to-one correspon- 
dence with the trace f G E^Ti.'^, and we have equivalences of norms 

WfW ^snp \\Ft\\^\\\tdtFt\\\. 

The corresponding reproducing formula Ft = e*'"^^'/, / = -FIr™ G E^lhi^Q is also valid 
for F solving the equations in R""^"*^, and the corresponding estimates holds. 

Proof, (i) implies (ii): As in the proof of Lemma [2.341 from Lemma [2.281 it follows 
that Yimt^oFt = f and limt^oo Ft = 0, and also that djEt = (-iT^D^e"*!^-^!/. 
Therefore Ft G C^{R+;n'') for all j. For j = 1, we get = {dt + {TsDEt = 
{dt -\-TB)Ft, since Ft G E'^'H^^. As explained in Section YTM this is equivalent with 
the two equations dt^xF = and dl^{BF) = 0. 

(ii) implies (i): The two equations can be written {—fi*dt + d*)BF = and 
{fidt + d)F = 0. Applying fi* and /i respectively to these equations and using 
nilpotence, we obtain fi*d*BF = and fidF = 0. Therefore Ft G 7i^, and since 
is closed we also have / G Hb- 

Next we write Ft = F^ + F^^ , where F^ := E'^Ft G E^Ti.'^ and similarly / = 
f^ + f~ . We rewrite the equations satisfied by Ft as dtFt + TBFt = 0. Applying the 
projections E^ to this equation yields 

dtF+ + \TB\F+ = 0, 
dtFi -\Tb\F- = 0, 
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since TsFt^ = ±\TB\Ft^. Fix t > 0. Then it follows that e^*-")!^^!^- jg con- 
stant for s E {t,oo) and that e^*~*'"^^'F^"^ is constant for s G (0,t). Taking limits, 
this shows that = lims^t+ e^^-^^l'^^lF" = lim.^oo e^*"")!^^!^^ = and F^^ = 
lim.^t- e(^-*)l^«lF+ = lim.^o e('-*)l^^lF+ = e-*l^sl/+. Therefore f = E E^n% 
and Ft = e-*!^^!/. 

To prove the norm estimates we note that ||/|| = lim^^o By Proposition 12.321 
g-t|TB| g^j-g uniformly bounded and thus sup^^Q < ||/||. Furthermore, using 

Proposition [231] with tpiz) = ze'^'^ shows that ^ ||/||. □ 

Remark 2.50. In proving that (ii) implies (i), it suffices to assume that \\Ft\\ grows 
at most polynomially when t — >■ oo. Indeed, from the equation Ff = e~*^*~*)'-^^lF~ 
for s > t, it then follows that 

||T^F-|| = (s-t)-'^||((s-t)rB)'^e-("-*)l^^lF-|| <Cis-t)-''\\F,\\ ^0, 

when s — * oo. Since Tb is injective, this shows that Ff = and therefore that 
Ft = e-*l^sl/ E E+n% as before. 

We now proceed by showing how, given data g in (Tr-i?^a^), we can solve for the 
trace / = -F|Rn by using the boundary operators Eb and A''^. 

Lemma 2.51. Assume that Tb satisfies quadratic estimates, so that the Hardy pro- 
jections E^ are bounded by Proposition let a"^ E C be given jump parameters 
and define the associated spectral point A := {a~^ + a~)/{a^ — a^). Then 

A — EBkNBk : Ti-B — ^ 
is an isomorphism if and only if the transmission problem (Tr-B^a^) is well posed. 

Proof. If we identify the fc- vector fields a;) in the transmission problem (Tr- 

B^a^) with the boundary traces f"^ and write f = f^ + f~ using Lemma [2. 49 j then 
we see that the transmission problem is equivalent with the system of equations 

N-,{a+E+, - a-E-,)f = N-,g. 

Using E'^f. = i(/ ± EBk) and adding up the equations, we see that the system is 
equivalent with the equation 

(A - EBkNBk) f = EBkg. 

a'^ — a 

This proves the lemma. □ 

Next we consider k = 1 and the boundary value problems (Neu-A), (Reg-A) and 
(Neu-'"-y4). By Lemma [2.491 we have the following. 

(1) (Neu-y4) is well posed if and only if the restricted projection : E\1H} — >■ 
NjTi.^ is an isomorphism. 

(2) (Reg- A) is well posed if and only if the restricted projection : E\7{^ — >• 
N^li} is an isomorphism, or equivalently if and only if A^"*" : E^lR} —>■ N^Ti^ 
is an isomorphism. Note that both N\ and project along N~1H} . 

(3) (Neu"'"-y4) is well posed if and only if the restricted projection : E^H} —>■ 
N'lH} is an isomorphism. 
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Proposition 2.52. Assume that T4 satisfies quadratic estimates. Then (Reg-A) is 
well posed if and only if (Neu^-A*) is well posed. 

Proof. We need to show that if A^"^ : E\l-l} —>■ N^IH} is an isomorphism, then so is 
A^^ : E\tlH} — s> N~7{^. The proof uses two facts. First that we have adjoint oper- 
ators (Ea)' = —Ea* and N' = N according to Proposition 12.171 and Lemma 12.471 
As in Remark I2.19[ this shows that we have dual spaces {E^H^ , E^^H^) a and 
{N+'H},N+'H})a, and we see that 

{N^f^g)A={f.9)A={f.E2,g)A. 

for all / G E\ii} and g G N'^'H^. Therefore, the restricted projections A^"*" : 
E+V} N+n^ and E^. : N+H^ E^M^ are adjoint. 

Secondly, if Rf and i?^ are two pairs of complementary projections in a Hilbert 
space as in Definition 12. 131 then : R^Ti RiTC has a priori estimates, as in 
Remark 12.181 if and only if i?^ : RfTi R^'H has a priori estimates. Indeed, both 
statements are seen to be equivalent with that the subspaces R^Ti and R^'H are 
transversal, i.e. that the estimate II/1 + /2II ~ ||/i|| + II/2II holds for all /i G R^l-i and 
f2 G R^H. To see this, assume that ||-Rf/2|| > II/2II holds for all /2 G R^H. Then 
11/2 II ^ \\Riifi + /2)|| < II /i + /2II for all fi G Rfn, which proves transversality. 
Conversely, assume that RfH and R^H are transversal. Then f2 — Rif2 = Rf f2 ='■ 
A G R+H for all /2 G R+H. Therefore pfMI = II/2 - /ill ~ II/2II + ||/i|| > IIMI- 
The same argument can be used to show that transversality also holds if and only 
if 11^2" /ill > ll/ill holds for all /i G RfH. 

To prove the proposition, assume that A^"*" : E^'W — N^'W is an isomorphism. 
It follows that the adjoint operator i?^* : N^ii} E^^Ti.^ also is an isomorphism. 
Using the second fact above twice, shows that 

E^ : N-n^ — > E^n^ and A^" : E^M^ — > N'n^ 

have a priori estimates. As these are adjoint operators as well, both must in fact 
be isomorphisms. In particular we have shown that (Neu"*"-/!*) is well posed. The 
proof of the converse implication is similar. □ 

When perturbing A, it is preferable not to have operators defined on spaces like 
E^Ti, which varies with A. We have the following. 

Lemma 2.53. Assume that Ta satisfies quadratic estimates. 

(1) If I — EaNa '■ Ti.^ —>■ Ti.^ is an isomorphism, then the Neumann problem 
(Neu-A) is well posed. 

(2) If I + EaNa '.ii} ^ii} is an isomorphism, or if I + EaN : ^ is an 
isomorphism, then the regularity problem (Reg-A) is well posed. 

(3) If I — EaN : — > Ti.^ is an isomorphism, then the Neumann problem 
(New^-A) is well posed. 

Proof. Assume for example that I + EaNa is an isomorphism. We need to prove that 
N^ : E^n^ N\ii^ is an isomorphism. Note that ii N\f = g where / G EXV}, 
then 

9 = Kf = W + NA)f = \{Ea + NA)f = \Ea{I + EANA)f. 
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If g e N^n\ let / := 2{Ea + NaY^Q- Then it follows that Q = N^g = \{Ea + 
NA){E^f), since iVj(E^ + Na) = \{Ea + Na - I - NaEa) = {Ea + Na)E^, so 
E2f = and therefore / G E^n\ 

A similar calculation proves well posedness of the other boundary value problems. 

□ 

Remark 2.54. More generally, letting k = 1 and {a~^,a^) = (1,0), i.e. A = 1, in 
Lemma 12.511 we see that / — EaNa is an isomorphism if and only if the restricted 
projections 

NX : E^ii} ^ NXii^ and 
: E+'fe ^ N^ii} 

are isomorphisms. Similarly, if (a+,Q;^) = (0,1), i.e. A = —1 in Lemma 12.511 
we see that / + EaNa is an isomorphism if and only if the restricted projections 
NX : E^n^ N^V} and A^^ : E^li} N^H^ are isomorphisms. 

We next turn to the Dirichlet problem (Dir-A), where we aim to prove an ana- 
logue of Lemma [2.491 which characterises the solution Ut as a Poisson integral of the 
boundary trace u. As discussed in the introduction, we shall use (Neu^-A) to con- 
struct the solution Ut- Given Dirichlet data u e L2(R"'; C), we form ucq G N~'U}. It 
then follows from Lemma 12.121 that the vector field Ft solving (Neu-'"-y4) with data 
mcq, has a normal component U := Fq which satisfies the second order equation 
(11. ip . We now define the Poisson integral of u to be 

Vt{u) := {Ft, eo), when Ft = e"*!^-^!/ and (/, eo) = u. 

Lemma 2.55. Assume that Ta satisfies quadratic estimates and that the Neumann 
problem (New^-A) is well posed. Let u G L2(R"; C) and let (0, oo) 3 t Ut{x) = 
U{t,x) G L2(R"; C) be a family of functions. Then the following are equivalent. 

(i) Ut = Vtu for allt>0. 

(ii) Ut G C2(R+;L2(R";C)), Vt,xUt G Ci(R+; L2(R"; 0"+^)) and U satisfies 
the equation 

diYt,^A{x)Vt,xU{t,x) = 0, 
and we have L2 limits \imt^o+ Ut = u, limt_,oo f^t = and linij^oo ^t,xUt = 0. 

If this holds, then UtyVt^xUt G C-' (R_|_; L2(R")) for all j > 1. Furthermore, Ut is in 
one-to-one correspondence with the trace u, and we have equivalences of norms 

\\u\\ ^ sup \\Ut\\ ^ |||t9tf/t|||. 

t>0 

Proof, (i) implies (ii): Assume that Ut = (Fj,eo), where Ft = e"*'"^-*'/ and (/, eo) = 
u. As in the proof of Lemma 12.341 from Lemma 12.281 it follows that Ft f, 
and therefore Ut u, and also that dtFt = (— iTAD-'e"*'"^-*'/. Therefore Ft G 
C''(R+; L2(R")) for all j, and so does Ut. For j = 1, we see that F satisfies the 
Dirac type equation since ITaI-F* = TAFt, and Lemma 12.121 thus shows that Ut 
satisfies the second order equation. Furthermore, we note from the expression (11. 5p 
for Ta, that V^Ut = {dtFt)^^. Thus V^f/j G C^(R+; L2(R"; C)) for all j, and yet 
another application of Lemma [2.281 shows that Ut = o(l) and Vt,xUt = o{l/t) when 
t —>■ 00. 
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(ii) implies (i): Assume Ut has the stated properties and boundary trace u. Con- 
sider the family of vector fields Gt '■= Vt^xUt- Since these satisfy Lemma F2.49( ii) for 
t > s > with boundary trace Gs, we obtain that Gs+t = e'^^'^^^Gg for all s,t > 0. 
For the normal components, this means that 

doUs+t = Vt{doUs), 

or equivalently that dsiUs+t — 'Pt{Us)) = 0. Since limg^oo f^s = 0, we must have 
Ug+t = VtiUg) for all s,t > 0. Letting s ^ 0, we conclude that Ut = Vtiu). 

The equivalence of norms ~ sup^^g ll^t II follows from the uniform boundedness 
of the operators Vt- For the equivalence ~ |||t5tt/t||| we use that (Neu-'--A) is 
well posed and the corresponding square function estimate for Ft from Lemma 12.491 
to get ^ 11/11 ^ \\\tdtFt\\\ ^ |||t<9t?7t|||, since for alU > we have \\dtUt\\ = 
\\N-{dtFt)\\ ^ WdtFtl □ 

We end this section with the proof of the non-tangential estimate ||A^^,(F)|| ^ ||/|| 
in Theorem ll.il 

Proposition 2.56. Assume thatTA satisfies quadratic estimates. Let Ft = e"*'"^-^'/, 
where f e E\7{^ . Then ||/|| ~ ||A^*(F)||, where the non-tangential maximal function 
is 

iV,(F)(x):=supf// \F{s,y)\^dsd'^^ , 

t>0 V JD{t,x) / 

and D{t,x) := {{s,y) E R"^^ ; |s — t| < Cot, \y — x\ < Cit}, for given constants 
Co G (0, 1) and Ci > 0. 

The proof uses the following lemma. 

Lemma 2.57. Let f E IH} and define Ht = {1 + itTsY^f G V} . Write Ht = 

(i) The normal component H}'^ satisfies the second order divergence form equa- 
tion 



[1 itdiv] 







" 1 " 


_a||x a||||_ 




itV 



hI'° = [1 itdiv] 



1,0 



flxx/ 



where we identify normal vectors ucq with scalars u, and the tangential com- 
ponent ifj^'" satisfies 



(ii) There exists p < 2 and q > 2 such that for any fixed tq < oo we have 



\fB{x,rot) J \fB{x,rot) J 

for all x e R" and t > 0. Eere M{f ){x) := supr>oiB(a; r) \fiy)\^y denotes 
the Hardy-Littlewood maximal function. 

Proof, (i) Multiplying the equation (1 + itTB)Ht = f with Mb we get 
(m{ix + ltd) + B-^{m{~^i* + itd*))B^Ht = Msf. 
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Similar to the proof of Lemma 12.121 we can now use the anticommutation relations 
from Lemma 12.71 to rewrite this equation as 

- (/i + itd)m + B-\-{-fj.* + itd*)m)B^Ht = Msf. 

since {m, fi + itd} = I, {m, —fi* + itd*} = —I and I — B^^IB = 0. Then apply 
(-/i* + itd*)B to obtain 

(/i* - itd*)B{fi + itd){mHt) = (-/i* + itd*){BN+ - N'B)/, 

since —fi* + itd* is nilpotent. Evaluating the scalar part of this equation, we get the 
desired identity. 

To find the identity for hI'\ we use the expression for Ta = ^bIt^i from Defini- 
tion Multiplying the equation (/ + itfA)Ht = / by AN+ - N'A yields 

a||||ifi^'" - aoo-f^/'°eo - it{AVHl'^ + d*fiAHt) = aim/^'" - aoo//'°eo. 

Evaluating the tangential part of this equation gives the desired identity, 
(ii) By rescaling, we see from (i) that it suffices to show that 



(2.13) 



\u{y)\'^dy 



1/9 



'B{x,ro) / \J B{x,ro) 

for all X G R" and all u and g satisfying an equation 



i/p 



\Vu{y)fdy] <M{\grf^{x), 



[l idiv] 



'o' a' 




' 1 " 




.^W^ ^\\\\. 




iV 


u 



\l zdiv] 



where A' is a matrix with same norm and accretivity constant as for A. Indeed, 
by rescaling we see that u{x) = g{x) = f{tx) and A'{x) = A{tx) satisfies 

this hypothesis. To prove fl2.13p . we use that the maps g ^ u and g ^ Vu have 
Lp(R'^) Lg(R") and Lp(R") — > Lp(R") off-diagonal bounds respectively, with 
exponent M for any M > 0, i.e. there exists Cm < oo such that 



(2.14) 



\u\\Lg{E) + ||Vm| 



<CM(dist iE,F)y\\g 



whenever E,F C R" and suppt? C F. To see this, let L := [l idiv] A' [l iV]*. 
Note that L : W^(W) Wa"^!^") is an isomorphism and that L : W^iR"") 
W~^(Tl"') is bounded. Then by the stability result of Sneiberg pi], it follows that 
there exists e > such that L : WplIV^) Wp^{TV^) is an isomorphism when 
|p — 2| < e. We then fix po E {2 — e, 2) and use Sobolev's embedding theorem to see 
that 



ll'^ll<?0 



+ I|Vm||po ^ hWw}^^ < \\Lu\\w~^i < \\g\\po. 



By choosing po close to 2, we may assume that go > 2. Thus we have bounded 
maps g u : Lp^^{F) — > Lg^^E) and g i-^ Vu : Lp^{F) — > Lp^{E), with norms 
< C. Also, by Proposition 12.251 the norms of g u : L2{F) L2{E) and 
g \—>- Vu : L2{F) —> L2{E) are bounded by (dist (i?, F))^^'^". Interpolation now 
proves (I2.14p for some Pq < p < 2, 2 < q < qq and M = Mq{1 — po/p)/(l — Po/2). 
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Finally we show how fl2.14p implies fl2.13p . Let E = Fq := B^x^Tq) and for k>l 
let Fk := B{x, 2Vo) \ B{x, 2'=-Vo). This gives 



/ HyWdy) + / \VuiyWdy) <J2^'''' ( \g{yWdy) 

JB{x,ro) J \JB{x,ro) J ;,^o ^"^^^ ^ 

oo ^ „ X 1/p 

<J22(n/p-M)kll \g{y)\Pdy] < (M( 1^?^) (x)) V^', 



provided we chose M > n/p. □ 

Proof of Proposition \2. 56\ To prove that ||A''^,(F)|| > ||/||, we calculate 

|iV,(F)f >sup / / / 

t>0 J R^J \y-x\<citJ \s~t\<cot 

sup-/ > sup ||F(i+co)tr ~ " ^"^ 

t>0 J|s_t|<cot t>0 



til ~ 



We have here used that F(^i^co)t = 6^^''^'^''°^*'^''^'^'^^ Fg, which by Proposition 12.321 shows 
that ||-F(i+co)t|| < ll-^sll when \s - 1\ < cot. 

To prove ||A^*(F)|| < ||/||, we note that since curli^i,.F = 0, we can write F = Vt,xU 
for some scalar potential f/, and we see that U solves the second order equation 
(II -ip . With notation D{t,x) = {{s,y) ; |s — t| < dot, \y — x\ < Cit}, for constants 



Co < Co < 1 and ci < ci < oo, and U '■=fjf)(ix^ U{s,y)ds dy, we have 
(II \tVUis,y)\'dsdy] ' < U-f \U{s,y) -U\^ds dy\ ' 

V JD(t,x) J V JDit^x) J 



< 



i/p 

\tVU{s,y)\Pdsdy " 

D{t,x) 



with 2(n+l) / (n+3) < p < 2. The first estimate uses Caccioppoli's inequality and the 
second estimate uses Poincare's inequality. Thus it suffices to bound the L2 norm of 
supj>o {§§D(t,x) W{s,y)Y' dsdyf'"^ . To this end, we write Fs = Hs + ipsiTA)f, where 
Hg := {I + isTa)'^ f and il){z) := e"'^' — (1 + iz)^^ . Using the quadratic estimates, 
the second term has estimates 



/ sup / \i)s{TA)f{y)Y dsdy] dx 

JYiP t>0 V JD(t,x) / 



< 



supf/ / \MTA)f{y)\'^]dx 

R" t>0 \J \s-t\<cot J \y-x\<cit 



Jo JR" J\y-x\<cis/{l-CQ) 



< rmTA)fr-<\\f\\' 

Jo s 



42 



PASCAL AUSCHER, ANDREAS AXELSSON, AND STEVE HOFMANN 



For the first term, we use Lemma I2.57( ii) and obtain 
' sup I 4 ■/ \Hs{y)\^ ds dy \ dx 

Rn t>0 V JD{t,x) J 

< [ sup sup f-f \Hs{y)\''dy\ dx 

7r" t>0 \s-t\<cot \Ib(x,cis/(1-co)) / 



<iiM(i/nii^/:<iii/ni^/:-"-"^ 



l2/p ~ II M I U2/p ll-l 112' 

using the boundedness of the Hardy-Littlewood maximal function on L2/p(R"). This 
completes the proof. □ 

3. INVERTIBILITY of unperturbed OPERATORS 

In this section we prove Theorem 11.41 and Theorem 11.11 for the unperturbed prob- 
lem, i.e. for Bq and Aq respectively. We do this by verifying the hypothesis in 
Lemma 12.511 and Lemma 12.531 i.e. we prove that T^k satisfies quadratic estimates 
and A — E^kN^k is an isomorphism, and that T^q satisfies quadratic estimates and 
/ ± EaqNao are isomorphisms respectively. 

3.1. Block coefRcients. In this section we assume that B = Bq E Loo(R"; >C(A)) 
have properties as in Definition 12.81 with the extra property that it is a block matrix, 
i.e. 

r^^x 

B, 



B 



in the splitting Ti = N Ti ® N^Ti. Note that B being of this form is equivalent 
with the commutation relations N^B = BN^. 

Lemma 3.1. Let B be a block matrix as above. Then 

Tb = t + b-^t*b, 

where T = Nmd = —iNd is a nilpotent first order, homogeneous partial differential 
operator with constant coefficients. 

Proof. Since N^B = BN^, it follows that Mb = N+ - B-^N-B = N+ - N- = N 
and 

Tb = N{md + B-^md*B) = T + B-^T*B, 

since A^^ = / and NB^^ = B~^N. The operator F is nilpotent since F^ = 
Nmd Nmd = NmNdmd = -NmNmd^ = 0. □ 

Remark 3.2. Note that if = F + F^, where T% = B-^T*B and F is nilpotent, is 
an operator of the form considered in [S], then 11^ intertwines F and F^ in the sense 
that UbTu = F^jH^M for all u G D(T*j^Ub) and n^F^jW = TUbu for all u e D(Fnij). 
Thus commutes with both F and F^ on appropriate domains. In particular, if 
pB = (1 + t^Ul)-^ and gf = tUB{l + t^n|)-i, then we find that TP^^u = P^Tu, 
T%Qfu = QfTu for all u e D(T) and F|jPf m = Pt^T*sU, FQf m = Qi^T%u for all 

u e D{r%). 

Theorem 3.3. Let B be a block matrix as above. Then 
(i) Tb satisfies quadratic estimates, and 
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(ii) EbNb + NbEb = and Nb = N. In particular X — EbNb is an isomorphism 
whenever A ^ {i, —i} with 

(A - EbNb)'' = T^(A - NbEb). 
A^ + 1 

Proof. For operators of the form F + B^^r*B, quadratic estimates were proved in 
[9] , with essentially the same methods as we use here in section 14. 1[ 

To prove (ii), note that since 5 is a block matrix, it follows that B^^N^Ti = N^Ti. 
Thus the projections associated with the splitting Ti = N^H © B^^N^H are 
A^^ and the associated reflection operator is A^^ = N. 

To prove invertibility of A — EbNb, we note that 

TbNb = (F + B-^T*B)N = -N{T + B-^T*B) = -NbTb, 

since N commutes with d, d* and B, and anticommutes with m. Thus 

EbNb = sgn{TB)NB = NBSgn{NBTBNB) = NBSgn{-TB) = -NbEb, 

since sgn(2;) is odd. Using this anticommutation formula we obtain 

(A - NBEB)i\ - EbNb) = A' + 1 - \{EbNb + NbEb) = A' + 1, 

and similarly (A — EbNb){X — NbEb) = A^ + 1, from which the stated formula for 
the inverse follows. □ 

3.2. Constant coefficients. We here collect results in the case when B{x) = -B G 
C{A) is a constant accretive matrix. In this case we make use of the Fourier transform 

acting componentwise. If we let 

then Tb, conjugated with JF, is the multiplication operator 

MJiO ■■= M-\zmfi^ - zB-'mf,lB)f{0, C e R". 
Lemma 3.4. For all t E H and ^ E R*^ we have 

\{tt + M^)-'\^it^ + \^\Y'^^- 

Proof. Let u = {it + M^)f. It suffices to prove that \\u\\^ > {f + I^HH/f • With 
B := mBm, it follows from the definition of that 

tuMbu = (F + B-^T*B)f, 

where F = i{t^ + ^^). Using Lemma [2.71 we get 

||(F + ngf = WTgf + ||r*^?ir = ((FT + T*T)g,g) = (t' + mUgf. 

Therefore our estimate follows from Lemma [2.21( ii). □ 

Proposition 3.5. // B{x) = B E C{A) is a constant, accretive matrix, then Tb 
satisfies quadratic estimates. 
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Proof. Using the lemma, we obtain the estimate 



1 + t2M| 



<t\M^\\{i-Mt^)-'\\{t + Mt 



~ l + t2|^|2- 



Thus using Plancherel's formula we obtain 



B 



1 + t^n 



u 



B 



dt 



where the last step follows from a change of variables s 




- m)vdi^\\u\ 



□ 



Next we prove that the Neumann and regularity problems are well posed in the 
case of a complex, constant, accretive matrix 



A 



)I and the action of T4 = Tb 



aoo Ooii 

Olio Ctiiii 

For this it suffices to consider B = I®A®I®I(£ 
on the invariant subspace li}. Recall that / G li} means that / is a vector field 
/ : R" ^ = C""*"^ such that (i/n = 0. On the Fourier transform side / G Ti^ is 
seen to correspond to a vector field / : R" such that ^ a /y = 0, i.e. / is such 

that its tangential part /y is a radial vector field. Thus the space 

H-n') = {/ G L2(R"; C"+^) ; e A /„ = 0} 

can be identified with L2(R"; C^) if we use {cq, ^/|^|} as basis for the two dimensional 
space to which /(^) belongs. Furthermore, the operator Ta is conjugated to the 
multiplication operator 



under the Fourier transform, where M4 := A^"*" — A A. We see from the ex- 
pression (11.51) for Ta that, at a fixed point ^ G S*""^, the matrix for in the basis 
{eo,0 is 



«00 



(aoii+ 0,10,0 



-I 

Theorem 3.6. Let A{x) = A E C{/\) he a constant, complex, accretive matrix. 
Then 

I ± EaNa -.n^ — ^n^ 

are isomorphisms. In particular, the Neumann and regularity problems (Neu-A) and 
(Reg-A) are well posed. 

Proof. From Remark 12.541 we see that it suffices to prove that all four restricted 
projections 

: E^n^ — > N^n^ 

are isomorphisms, or equivalently that the constant multiplication operators 
are isomorphisms on T{1H}). For the projections x±{^i) conjugated to the Hardy 
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projection operators E^, we observe that Mt^ = tM^ and hence X±(^^c) = X±iM^) 
for alH > 0. Thus it suffices to verify that 



are isomorphisms for each ^ E S"' ^. For such fixed ^, using the basis {eo,^} from 
above, N^Hl = {/ e Hi ; cq a / = 0} is spanned by [l O]* and N+Hl = {/ G 
; (A/, Co) = 0} is spanned by [(aoi|,0 ~Q-oo]*- Indeed {A{zeo + w^) , Cq) = 
2;aoo + w(aoi|,0- 

If we call and the two eigenvectors of M^, it follows that X±(^c)'^^ ^^re 
spanned by these, as these subspaces are one dimensional. It suffices to show that 
[l O] and [(aoii,0 ~'^oo] are not eigenvectors of M^. We have 



^(«0|| + 0,10,0 



(«0||,0 

—aoo 



1 



(«oii + a||0,0(«oii,0 - 

-(«0||,0 



Clearly the normal vector is not an eigenvector. To prove that the second is not an 
eigenvector, note that the cross product is 

-(«oii,0^ + («oii + a||o,0(«oii,0 - «oo(a|iii^,0 = («iio, O(«oii, - «oo(a|iii^, 0- 
The right hand side is non-zero since 



0,10 a. 





z 




z 






5 





(aiio,0 (aiiii^.O 





z 




z 




w 




w 



is a non-degenerate quadratic form as A is accretive. 



□ 



Remark 3.7. We note that the method above also can be used to show that / ± 
EaN : Ti.^, with the unperturbed operator N, are isomorphisms when A is 

constant. Here we also need to observe that the tangential vector [O l] is not an 
eigenvector to M^. 

3.3. Real symmetric coefRcients. In this section, we assume that B* = B. We 
ffist prove a Rellich type estimate. 

Proposition 3.8. Assume that B* = B and that f G E'^Vb or f e E]^Vb- Then 

{Bf, /) = 2 Re (eo j (Bf), eo j /) = 2 Re (eo a (Bf), eo a /) . 
In particular \\f\\ ^ \\N^f\\ ^ \\N^f\\ ^ ||iV+/|| ^ ||iV+/||. 



Proof. It suffices to consider / G E^Vb as the case / G E]^Vb is treated similarly. 
We use Lemma 12.341 and write Ft := e"*l^slj. Hence (0, oo) 9 t ^^ G is 
different iable, limj^o-^* = / and lim^^oo -^t = 0. Furthermore Ft G D{dx) and 
BFt G D(d*) for all t G (0, oo). We note the formulae 



dt,x'>T^Gt + mdl^Gt 
Bmdt.T.Ft 



t,x-rt 



-dtGt, 
-md*t BFt. 
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The first identity, which we apply with Gt = BFt, follows from Lemma [213 whereas 
the second is equivalent to dfFt + TsFf = 0, and follows from Lemma .341 We get 

{Bf, f) = - {dtBFt, Ft) + {BFt, dtF^) = -2 Re / {dtBFt, F^) 
Jo Jo 

POO POO 

= 2 Re / {dl^mBFt + mdl^BFt, Ft) = 2 Re / {dl^mBFt - Bmdt,^Ft, Ft) 
Jo ' ' Jo ' 

POO 

= 2 Re / {dl^mBFt, Ft) - (mBFt, dt,^Ft) 
Jo 

POO 

= 2 Re / {{d* - fi*dt)mBFt, Ft) - (mBFt, (d + fidt)Ft) 
Jo 



= -2 Re / {dtmBFt, fiFt) + (mBFt, dtfiFt) 
Jo 

= 2 ReimBf, fif) = 2 Re(eo a (Bf), eo a /). 

Note that all integrals are convergent since < min(l,t~'') and since dtFt = 

-TeFt with HT^FJ < min(t^-i,r^) if / e ^[\TbY) n R(|Tb|-^). This follows as in 
the proof of Lemma F2.34[ Furthermore, we note that 

(eo A Bf, eo A /) = {Bf, cq j (eo a /)) 

= {Bf, / - eo A (eo J /)) = {Bf, f) - (eo j {Bf), e^^f). 

Together with the calculation above this proves that {Bf, /) = 2 Re [eoj{Bf), eoj/). 

To prove that ||/|| ^ ||iV+/|| ^ ||iV^/||, it suffices to show that ||/|| < ||iV|/|| 
since are bounded. From the Rellich type identities above we have ^ 
||eo J (i?/)||||eo J /II < IliV^/ll 11/11 which proves ||/|| < ||iV^/||, and using the other 
identity we obtain H/f < ||eo a (5/) || ||eo a /|| < ||/|| ||iV+/||. The proof of ||/|| ^ 
11^^/11 - WNbIW is similar. □ 

We note that Proposition 13.81 also proves that the norms of the two components 
of / G E^Vb in the unperturbed splitting H = N'H © N'^H are comparable. 

Corollary 3.9. Assume that B* = B and that f G E^Vb or f e E^Vb, and 
decompose B in the splitting Ti = N^Ti. © N'^H as 



B 



B±± i?x 
Bu I Bn 



Iffii=N+f andf, = N-f, then 

{B^J\,f^) = {B,,f„f,). 
In particular \\f\\ ^ \\N+f\\ ^ \\N'f\\. 

Proof We obtain from Proposition ES] that Re{B^J^ + B^J^^,f^) = Re{B^^J^ + 
i?iiii/ii, /||). The corollary now follows since -Bj_x, -Bim > and = -B||±. □ 

Next we turn to quadratic estimates for the operator Tb- As before we assume 
that B is as in Definition 12.81 and that B = B*. For the rest of this section we shall 
also assume that 

(3.1) Bi = Bi,=0, for J > 2. 
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Note that in particular this is true if 5 = / ©A©/©...©/, where A G 
Loo(R";i:(Ai)). 

To prove quadratic estimates, we shall use Proposition l2. 361 where we verify the hy- 
pothesis separately on the subspaces and Hb, which is possible by Lemma B.46[ 



Lemma 3.10. Assume that B* = B and h3. 1\) holds. Then 

11/11 < \\\tdtFt\\\^, for all f e E^Vb, 

where Ft = e^*'^^'/. 

Proof. As the two estimates are similar, we only show the estimate for / G E^Vb- 
Note that F{t,x) satisfies fl^TT^ . i.e. 

dF = cq J dtF, 
d*{BF) = -eoA{BdtF). 

Splitting both sides of both equations into normal and tangential parts, with nota- 
tion as in Corollary 13.91 we obtain 

(3.2) rfi^i = mdtF^, 

(3.3) dF^ = 0, 

(3.4) d*{B^^F^ + B^,F,) = -m{B,^dtF^ + B,,dtF,), 

(3.5) d*{B,^F^ + B,,F,) = 0. 

We have here used that eo-i/ = mf if / is normal, and cqa / = mf if / is tangential. 
A key observation is that the first term on the right hand side in (13. 4p vanishes since 
B^^^Fj_ = 0. To see this, note that 

= Bl^Fl + Bl^Fl + ... + B'^^'^F^^K 

By hypothesis (13. ip . B'^^j_ = . . . = Bl^^'^ = 0. Furthermore, writing F| = FqCq with 
the function Fq being scalar, we get from (13.30 that = d^FQCo) = (VFq) a Cq, so Fq 
is constant and therefore vanishes, and thus so does Equation (13.41) reduces to 

(3.6) ^iiii^t^i = ~md*{B^^F^ + ^^yF,,). 
We calculate 

POO 

'<(5|iii/ii,/ii) = -2Re / iB,,dtF„F,)dt 

Jo 

POO 

= 2Re/ [d*{B^^F^ + B^\\F^;),mF\^dt 
Jo 

= -2Re^ ((c/*(5^x9tF^ + 5^ii9ti^i),mF||) + {d*{B^^F^ + B^^^F^^),mdtF^^)y dt 
= 2Re^ (^{B^AF^ + B^,dtF,,dtF^) + {mB^^^^dtF^^,mdtF,)y dt 

POO 

= 2 Re / (BdtF, dtF)t dt ^ \\\tdtF\f. 
Jo 

Here in the third step we use (13.61) . in the fourth step we integrate by parts, in the 
fifth step we use duality and that dmF\\ = —mdF^^ = —dtF^ by (13.21) for the first 
term and again (13. 6p for the second term. Finally in step 6 we use that m is an 
isometry and that B»^Fj_ = 0. □ 
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Next we turn to the subspace ii} . 

Lemma 3.11. Let A G Loo(R"; /^(A"*^)) he real symmetric and A > k > 0, and let 
B = I®A®I®...®I. Then 

\\f\\<\\\tdtF\U, for all f e E%Vl, 

where F = e^f"*'"^^'/. 

Proof. Recall that if / G E^V}^ and F = FqCq + i^i = e^*''^^'/, then Lemma [2ll2] 
shows that Fq satisfies the equation 

divt,c,A{x)Vt,xFoit,x) = 0. 

Therefore, by the square function estimate of Dahlberg, Jerison and Kenig [13] and 
the estimates of harmonic measure of Jerison and Kenig [17], we have estimates 
ll/oll ^ |||^Vf,x-Fo|||±- Hence applying the Relhch estimates in Proposition 13.81 we 
obtain 

11/11 ^ ll/oll < pv,.Foiiu < mF\u, 

since 9jFo = doFi, as dt^xF = 0. □ 
We are now in position to prove the main result of this section. 

Theorem 3.12. Let A G Loo(R-"; '^^(A"'^)) be real symmetric and A > k > 0, and let 
B = I®A®I®...®L Then Tb has quadratic estimates, so that in particular 
Eb = sgn(TB) : Ti. —>■ Ti. is bounded. Furthermore we have isomorphisms 

I ± EbNb -.n — >n. 

Proof. To prove that Tb has quadratic estimates it suffices by Proposition 12.361 to 
prove 

\\f\\<\\\tdtF\U, ioi ^\\ f e E^Vb. 

To this end, we split 

/ = /O + /l + • • • + /n+l + / 

where fj G Vb and f gVb- Lemma [2.461 shows that e^*'"^^' preserves these subspaces 
so that similarly 

F = Fo + Fi + ... + Fn+i + F, 
where {Fj)t G Vb and Ft G Vb- It thus suffices to prove that < |||)f:(9iFj|||-|-, 

J = 0,l,...,n + l, and 11/11 < \\\tdtF\\\±. Lemma Emi shows that ||/|| < \\\tdtF\\\± 
and Lemma [3.111 shows that ||/i|| < |||t(9tFi|||±. Furthermore, for j 7^ 1 we observe 
that Tb = Tb = Tj on the subspace V^, where / denotes the identity matrix. Thus 
Fj = e^*l^^l/j- and it follows from Proposition 13.51 that \\fj\\ < \\\tdtFj\\\±. 

Applying Proposition 12.361 and Proposition 12.321 now shows that Eb = sgn(TB) : 
7i ^ 7i is a bounded operator. To show that / ± EbNb : 7i — > 7i is invertible, note 
that the adjoint with respect to the duality (-, ■)b from Definition 12. 161 is 

/ T NbEb = Nb{I T EbNb)Nb, 

according to Proposition 12.171 Having established the boundedness of Eb, the Rel- 
lich estimates clearly extends to a priori estimates for all eight restricted projections 

: E^n — > N^n, : E^H — > N^H. 

As in Remark 12.541 this translates to a priori estimates for / ± EbNb and their 
adjoints, which proves that they are isomorphisms as in Remark 12. 181 □ 
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Remark 3.13. Using instead Corollary I3.9[ we also prove as in Theorem 13.121 that 
/ ± EbN : Ti — > Ti, where N is the unperturbed operator, are isomorphisms. 



4. Quadratic estimates for perturbed operators 
In Section [3] we proved that Tbq satisfies quadratic estimates 

(4.1) IIIQfVIII ^ 11/11, fen, 

for certain unperturbed coefficients Bq: 
(b) Block coefficients 



Bo 



(So)xx 

(Bo) 



(c) Constant coefficients Bo{x) = 5o, x G R", of the form Bq = /©Aq©/©. . .©/, 

i.e. Bq only acts nontrivially on the vector part, 
(s) Real symmetric coefficients of the form Bo = lQ)AoQ)I®---Q)I- 

Note that for (c), we did prove quadratic estimates for general constant coefficients 
Bq in Proposition 13. 5[ However, since we only prove invertibility I ± EbqNbq on the 
subspace H.^ in Theorem 13.61 we shall only prove perturbation results for constant 
coefficients of the form = / © ©/©...© /. 

In this section we let Bq G Loo(R"; >C(A)) be a fixed accretive coefficient matrix 
with properties (b), (c) or (s). Constants C in estimates or implicit in notation < 
and ~ in this section will be allowed to depend only on ||5o||oo, f^Bo dimension 
n. Note that this is indeed the case for the constants implicit in fl4.ip . 

We now consider a small perturbation B G Loo(R"; '^^(A)) of Bq: 

\\B — -Bo||loo(R";£(a)) < eo- 

Throughout this section, we assume in particular that eo is chosen small enough so 
that B has properties as in Definition 12.81 with ||-B||oo ^ 2||i?o||oo and > ^kbq- 
The goal is to show that 

(4.2) lligf /III < C\\f\\, whenever \\B - Bq\\^ <e, feH, 

where C = C(||i?o||oo, /^Bq' '^) ^ — ^{\\Bq\\oo, i^Bq:^) ^ ^o- Since the properties 
of Bq are stable under taking adjoints, we see that the quadratic estimate (12. 7p 
follows from (14. 2 p and Lemma 12.351 In order to prove (14. 2 p we make use of the 
following identity, where we note that each of the the six terms consists of three 
factors, the first being one of the operators from Corollary 12.231 the second being a 
multiplication operator £ with norm ||£^||oo ^ and the last factor being one of the 
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)Bo 

t 



operators from Corollary 12.231 but with B replaced by Bq. 

(4.3) gf - gf " = pf ((M^i - m^^)Mb,) gf ° 

(4.4) - zPf [m^\B-' - B,')Bo) {td%Pf" 

(4.5) -t[pfM^Hd%) [b-\B - B,))PI 

(4.6) - gf ((M^i - M^;)Mbo) {tTs^Ql 

(4.7) + zgf (M^i(fi-i - B,^)B,) (tdr^^Ql 

(4.8) + z(gf M^4c?^) (5-1(5 - i?o))gf' 
Recall that 

Tb = -iM^\d + B-^d*B), where d = imd, Mb = - B'^ B , 

pB = (1 + t^Tl)-! and Qf = tTB^l + t^Tl)-^ and similarly for Eq. The identity 
(14.3114. 81) is established by using 

F(l + _ x(l + X2)-i = (1 + Y^)-\{Y -X)- YiY - X)X)(1 + X^)-^ 

with y = ITb and X = tTe^, and then inserting 

Y-X = Tb-Tb, = (M^i - M^l)MB,tTB, 

- i{M^\B-' - B,')Bom%^) - i{M^H<r^){B-\B - B,)). 

Lemma 4.1. Assume that for all £ G Loo(R"; '^^(A)) and all operator families 
{Qt)t>o with L2 off-diagonal hounds \\Qt\\off < C as in Definition 2.24, we have 
the two estimates 

(4.9) Ill4^(trfgf')lllop<c(|||g,lllop + i)ll^lloo, 

(4.10) \\\QtmBoQ?')\\\ov < ciwmwop + i)ii^iioo. 

Then there exists e > such that Ili4.2\ ) holds. 

Proof. First consider the terms (14. 3p . (14. 4p and (14.81) . Using the uniform boundedness 
of Pf and QfM~^td*B from Corollary E^Sl we deduce that 

IIKSDIIlop < ||5-5o||oo|||gfl||op, 

|||gaD|||op<||i?-5o||oo|||tC^BoPf°|||op, 

|||(!MD|||op<||fi-5o||oo|||gfl||op. 

Observe that td^^Pf'' = tPl^MBoQf° as in Corollary ESSl where the Hodge pro- 
jection P^^ and Mbq are bounded. Thus we see from (14. ip that 

MWop + IIKSaDlllop + llldMDIIlop < \\B - 5o||oo. 



To handle the terms (l4.5H4.7p we introduce the truncated operator families Q] : = 
xif)Qt S'lid Ql := x(t)Pt^M^^td*B, where xit) denotes the characteristic function of 
the interval [t~^, t] for some large r. Note that Pf M^^td*^ = iQfP^ as in Corol- 
lary [2231 and therefore IHg^ |||op = IligtPslllop ^ llig* lllop. To use the hypothesis on 
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the terms fl4.5ti4.7l) we note that the last factors are 

tTs^Qf" = -zM^^itdQf") - tM^^{td%^Qf°), 

-^tdrJ,^Qf' = -Ktdi,Qf') 

respectively. Thus we get from (14. 3114 .8 1) , after multiplication with xif)^ that 

\\\xit){Qf - Qf")lllop < C(|||x(t)Qf lllop + i)l|5 - ^olloo, 

and thus if \\B - Bq\\^ < e := 1/(2C) that 

iiixQf iilop < c\\B -B4^ + mpu < c\ 

since |||x(t)Qf |||op < c)o. Since this estimate is independent of r, it follows that 

\\\Q?f\\\ < 11/11- □ 

Before turning to the proofs of fl4.9p and fl4.10p . we summarise fundamental tech- 
niques from harmonic analysis that we shall need. We use the following dyadic 
decomposition of R". Let A = Uj°l-oo where := {2^{k + (0, 1]") : G Z"} 
if 2^^^ <t< 2K For a dyadic cube Q G A23, denote by 1{Q) = 2^ its sidelength, by 
\Q\ = 2"-' its volume and by Rq := Q x {0,2^] C R"^^ the associated Carleson box. 
Let the dyadic averaging operator At : Ti. ^ Ti he given by 

Atu{x) := Uq := / u{y) dy = [ u{y) dy 
Jq IvI Jq 

for every x G R" and t > 0, where Q G A^ is the unique dyadic cube that contains 

X. 

We now survey known results for a family of operators 0^ : 7i ^ 7i, t > 0. For 
the proofs we refer to [3J and pj. 

Definition 4.2. By the principal part of (Ot)t>o we mean the multiplication oper- 
ators 7j defined by 

jtix)w := {Qtw){x) 

for every w G A. We view w on the right-hand side of the above equation as the 
constant function defined on R" by w{x) := w. We identify 7t(x) with the (possibly 
unbounded) multiplication operator '■ f{x) t— > 'yt{x)f{x). 

Lemma 4.3. Assume that Qt has L2 off-diagonal bounds with exponent M > n. 
Then Qt extends to a bounded operator L^o ~^ Ll^'^. In particular we have well 
defined functions jt ^ Li^^CR^; C{A)) with bounds 

\ltiy)\' dy < WQtWlff 

for all Q E At. Moreover ||7tA|| ^ l|0t||oif uniformly for all t > 0. 

We have the following principal part approximation Qt ^ It- 

Lemma 4.4. Assume that Qt has L2 off-diagonal bounds with exponent M > 3n 
and let Ft : R" ^ A be a family of functions. Then 

jtAt)Ft\\\<\\Qt\\offmFt\\\, 
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where V/ = V ® / = X]j=i ® i^jf) denotes the full differential of f . Moreover, if 
Pt is a standard Fourier mollifier (we shall use Pj = (1 + t'^Il'^)~^ where 11 = F + F* 
and F is an exact nilpotent, homogeneous first order partial differential operator with 
constant coefficients as in [9] j and Ft = Ptf for some f E H, then |||7iy4((P(— /)/||| < 
l|0t||oif||/|| and \\\tVPt\\\op<C. Thus 

\mPt-itAt)f\\\<\moff\\fi 

Definition 4.5. A function 'j{t,x) : R"^"*^ ^ A is called a Carleson function if there 
exists C < oo such that 

|7(t,x)p^<C7^|g| 

for all cubes Q C R". Here Rq := Q x (0,/(Q)] is the Carleson box over Q. We 
define the Carleson norm ||7f||c to be the smallest constant C. 

We use Carleson's lemma in the following form. 

Lemma 4.6. Let'jtix) = 'j{t,x) : H^'^ A be a Carleson function and let Ft{x) = 
F(t,x) : R"^""^ ^ A be a family of functions. Then 

WbtFtlW < ||7t||c;||iV*(i^*)ll, 
where N^,{Ft){x) := sup|^_3,|^jFt(?/)| denotes the non-tangential maximal function 
of Ft. In particular, if Ft = Af for some f e H, then ||iV,(A/)|| < ||M(/)|| < 
11/11, where M{f ){x) := sup,,>oiB(^ \ f{y)\dy denotes the Hardy-Littlewood maximal 
function, and thus 

IhtA/lll < Il7*llcll/I|. 

Lemma 4.7. Assume that 0^ has L2 off-diagonal bounds with exponent M > n. 
Then 

\\Qtf\\c<{\\\et\\\op + \moff)\\f Woo, 

for every f G Loo(R"'; A). In particular, choosing f = w = constant we obtain 

WltWc ^ ll|0t|||op + l|0t||o#- 

4.1. Perturbation of block coefficients. In this section we assume that Bq is a 
block matrix, i.e. we assume that 



'(5o)xx 

(Po)ii 

in the splitting Ti. = N~T-C © N^Ti. Our goal is to prove Theorem 11.41 by verifying 
the hypothesis of Lemma 14.11 We recall from Lemma 13.11 that 

= n^o = F + Bq^V*Bq, where F := -iNd, 

is an operator of the form treated in [9]. Thus with a slight change of notation for 
we need to prove the following. 

Theorem 4.8. If Bq is a block matrix and if {Qt)t>o is an operator family with L2 
off-diagonal bounds \\Qt\\off^ C as in Definition \2.24\ then with notation as above 

(4.11) \\\QtS{t'rr^X')\\\op < ciWmWo, + i)\\s\U 

(4.12) |||4^(t^F^/Pf°)|||„, < C(|||4|||o, + l)||^||oo, 
where P^° = {l + t^Ul^)~\ 
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We note that if (14.111) holds with T replaced by T* and with Bq replaced by -Bq ^, 
then also (14.121) holds. This follows from the conjugation formula 

Thus it suffices to prove (14. lip , as long as we only use properties of (T,Bq) shared 
with (T*,Bq^). To this end, we let 0^ be the operator 

and denote by 7t(x) its principal part as in Definition 14.21 We note that we have a 
Hodge type splitting H = N(T) ® U(T%^) by Lemma [2l2l](i), and since 6t|N(r|jj = 
it suffices to bound |||6(/||| for / G N(r). We do this by writing 

(4.13) Qtf = Qt{I - Pt)f + {QtPt - ltAt)f + jtAf, 

where 11 := F + F* is the corresponding unperturbed operator and Pt := {l + t'^lP)~^ 
and Qt :=m(l + t2n2)-i. 

Lemma 4.9. We have, for all f G A/(F), the estimate 

|||e,(/-p,)/|||<c(|||4lllop + i)ll^llooll/l|. 

Proof, life N(F), then {I-Pt)f = tVQtf G N(F), which shows that t'^TT*B^Pf° {I - 
Pt)f = (/ — Pf"){I — Pt)f ■ To prove the estimate, we write 

- Pt)f = QMI - Pf'Wtf = Qt£f - QtSPtf - Qt£Q?'iJiBl'^)Qtf. 

where we recall that = tHs^Pf". Clearly < Ill4lll ll^l|oo||/||- For the 

second term, we write := Qt£ with principal part 7* as in Definition 14.21 and 
estimate 

|||e*pj|||<|||(0*p,-7*A)/||| + |||7*A/lll 

< ll^llooll/ll + ||7*l|c;||/|| < ll^llooll/ll + (llig^lllop + l)||^||oo||/||. 

In the second step, we used for the first term that ||6t||off ^ ll^lloo in Lemma 
and for the second term we used Lemma 14.61 In the last step, we used Lemma 14.71 
on the last term and that |||0t|||op ^ |||<5t |||op||i^||oo- Finally we note that 

|||4^gf"(n^;F)g,/||| < ||^||ooll|g*/lll < ll^llooll/ll, 

since the Hodge projection H^^^F is bounded by Lemma [2.21( i). □ 

To estimate the last term in (14.131) we shall apply a local T{h) theorem as in [21 E]- 
We here give an alternative construction of test functions to those used in [9], more 
in the spirit of the original proof of the Kato square root problem [3]. 

Lemma 4.10. Let T he a nilpotent operator in Ti, which is a homogeneous, first 
order partial differential operator with constant coefficients. Denote by Ar C A the 
image of the linear functions u : R" — > A under T , where we identify A with the 
constant functions R*^ — >• A. 
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Then for each w G Ar with \w\ = 1, each cube Q C R" and each e > 0, there 
exists a test function /^^, G H such that /^^ G R{T), \\f^J < \Q\^/^, 



II < _l_|g|i/2 



I rSQj (, 



el{Q)\ 



and 



JQ 



W 



< .1/2 



Proof. Let u{x) be a linear function such that w = Tu and supgg |u(a;)| < 1{Q), and 
define ifg := T{riQu), where rjg is a smooth cutofi^ such that riQ\2Q = 1, supp (r/g) C 
3Q and ||Vr/Q||oo < l//(<5)- It follows that 

Wq eR(r), wq\2q = w, suppwqCSQ and ||wq||oo < C*. 

Next we define the test function fg^^ := P^°wq, where we write / = 1{Q). Using 
Corollary it follows that ||/^J| < \Q\^/^ and ||r^J-J| < ^^\Q\"^ and since 
r commutes with -P^^f", it follows that /g ^ G R(r). To verify the accretivity property, 
we make use of [9], Lemma 5.6] which shows that 



JQ \JQ \JQ 



< .1/2 



1/4 



1/4 

So„„_|2^ < pl/2 



where we used that wq G R(r) in the second step. 



□ 



Proof of Theorem 4-S We have seen that it suffices to prove fl4.1ip . and to bound 



each term in fl4.13p for / G N(r). The first term is estimated by Lemma [4.91 and 
the second by Lemma 14. 4[ To prove that the last term has estimate |||7ty4j/||| < 
(IIIQtlllop + l)||^||oo||/||, it suffices by Lemma 1^16) to prove that 

ll7*l|c<(|||4ll|op + l)||^||oo. 



To this end, we apply the local T{b) argument and stopping time argument in P, 
Section 5.3]. Note that R(r) C L2(R'^; Ar) and thus, since F is an exact nilpotent 
operator, it follows that N(r) C L2(R'';Ar). Furthermore Atf G L2(R";Ar) if 
/ G L2(R";Ar). Thus it suffices to bound the Carleson norm of 7f(x) seen as a 
linear operator 7f(x) : Ar — A. The conical decomposition [jj^^yK^, of the space 
of matrices performed in P, Section 5.3], here decomposes the space C{Ar; A) and 
for the fixed unit matrix u G C{Ar; A) we choose w G Ar and w E A such that 
\w\ = \w\ = 1 and z/*(w) = w. With the stopping time argument in [9, Section 5.3], 
using the new test functions fg from Lemma I4.10[ we obtain 

htrc< sup ^ // Ux)AJ-{x)\'^, 

where fQ=fQ^ for a small enough but fixed e. To estimate the right hand side we 
use (14.131) with / replaced with the test function fn, estimate the first two terms 
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with Lemma [4.91 (which works since /g G R(r)) and Lemma F4.4[ and obtain 
'Rq 



<[[ Mr,-ej-y-^+ [[ |e./,n^^ 



'Rq JJRq 
Using that \\V*j,J^\\ < j^^W^ we then get 

II©*/qII = IIQ*^(irpf°)t(r* /-)|| < \\sUjh\Q\'" 



dxdt 



This yields 



Rq 



005 



which proves that //^^ |7tA/^P^ < (IIIQ* lllop + □ 

Remark 4.11. (i) Note that using the new test function from Lemma [4. 101 sim- 
phfies the estimate of the term {'ytAt — Qt)fQ in the above proof as compared 
with the proof of P, Proposition 5.9]. The useful new property of the test 
functions from Lemma [4. 101 is that they belong to R(r). 
(ii) In the proof of Theorem 14. 8[ we only estimate the Carleson norm of the 
restriction of the matrix 7^ to the subspace Ar as this suffices since we want 
to bound the quadratic norm of 'jtAtf, and A^f is always Ar valued. However, 
to prove (14.111) and (I4.12p . we use F being either Nmd or Nmd* . In these 
two cases, the space Ar is either the orthogonal complement of span{l,eo} 
or span{eo,i,...,n, ^i,...,n\ respectively. Note also that block matrices preserve 
these spaces Ap. It is seen that in these two cases 74 = on Ap, so actually 
we do get an estimate of the Carleson norm of the whole matrix 7^. 

4.2. Perturbation of vector coefficients. In this section we assume that the 
unperturbed coefficients Bq are of the form 

5o = / © Ao ©/©...© /, 

i.e. Bq only acts nontrivially on the vector part, and that Aq is a matrix such that 
has quadratic estimates. Note that this hypothesis is true if Aq is either real 
symmetric, constant or of block form, by Theorem 13. 12^ Proposition 13.51 and Theo- 
rem 13.31 respectively. Our goal is to prove Theorem II. 1^ by verifying the hypothesis 
of Lemma HH], as well as proving Theorem ll.3[ We start by reformulating Lemma ITT] 
in terms of e~*'"^^o', acting only on functions / in one of the two Hardy spaces E^^Ti,, 
instead of Pf° . 

Theorem 4.12. If Bq is as above and if \\Qt\\off < C, then for all f G E^^H we 
have estimates 

(4.14) \\\QMFt-f)\\\<C{\\mip + l)\\Sl 



(4.15) 



QtSd / Esds 







<ci\\m\\op+i)\\£\ 
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where Ft := e ^l^^olj ffi^ extension of f as in Lemma \2.4^ The corresponding 
estimates for f G E^^Ti also holds. 

Proof that Theorem \4.1S\ implies Ili4.9\ ) and Ili4.10\ ). We first note that it suffices to 
prove fl4.9l) and fl4.10p for all / G E^^^H and all / G E^^H since we have a Hardy 
space splitting H = E^^Ti. © E]^^T-C. We only consider / G E^^H since the proof for 
/ G E^^Ti is similar. 

Now let / G E^^H and use Proposition 12.311 which shows that if G "^{3°), then 
\\MTbo)\\\op < IIIQflllop < C. For the estimate gl]), we write 

QtStdQf^^ = Qt8{dT^l){I - e-*l^-ol) + Qt8{dT^l)MTB,)^ 

where i){z) = e'^^^ - (1 + z'^)'^. Note for the ffist term that T^l{I - e"*l^sol)/ = 
Tj^^ (f — Ft) = — TgJ^ dsFgds = F^rfs. Therefore Theorem H?T2l the boundedness 
of dT^^ and Proposition 12.311 give the estimate (14.91) . 
For the estimate fl4.10p . note that it suffices to estimate 

QtStTs^Qf' = QtS{I - Pt""") 
since iMbqTbq = d + d^^^. But this follows immediately from fl4.14p since 

Qt£{I - Pt''")f = Qt£{f - Ft) + Qt£MTBo)f, 
with the same as above. □ 

Remark 4.13. In the case when Bq is a constant matrix, we can estimate Qt£{I — 
Pt^°)f directly. We prove that 

|||Q,^(/-Pf°)/|||<(|||4ll|op + l)||^||oo||/|| 

as follows. Clearly |||Qt^|||op ^ |||Qt|||op||^||oo- For the second term we write Qt '■ = 
Qt£Pf°- Inserting a standard Fourier moUifier Pj, we write Bj = 0t(/ — Pt) + &tPt- 
Here \\\et{I - Pt)f\\\ < \\\Pf°{I - Pt)f\\\ < ||/|| is easily verified using the Fourier 
transform. On the other hand 



llie.pjiii < |||(e,p, - 7t^)/||| + III will < 11/11 + htWcl 

using Lemma 14.41 and Lemma 14.61 However, since Bq is constant, we have that 
TbqW = if w is a constant function, and therefore 

-ftw = Qt£{P^'w) = Qt£w. 

Therefore Lemma Wl\ shows that ||7t||c ^ (il|Qt|||op + l)||^||oo- 

We now set up some notation for the proof of Theorem 14.121 We decompose the 
function Ft := e"*l'^sol/, where f e Ft H, 



as 



(4.16) Ft = P° + (P/'" + P/'") + F^ + ... + F^+\ 

and similarly for / = limi^o+ Ft. It is important to note the special property that the 
normal component of the vector part P/'^ = P/'°eo has by Lemma 12.121 it satisfies 
the divergence form second order equation diV(^a.Ao(x)Vt,x'P^'*' = 0. Furthermore, 
we decompose the matrix Aq as 

Oj\\ I flllll 
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in the splitting Ti = N^H © N^H. We view the components aj_j_, a^y, Oyx and 
a|||| as operators, and write a^^{f^'^eo) = (aoo/^'°)eo, Oxn/^'" = (aoy ■ /^'")eo and 
a||x(/^'°eo) = /^'"o-iio, where scalar and aoy and ayo are vectors. 

We introduce an auxiliary block matrix 



I® An 







Lemma 4.14. Let Ft := e-'l^^ol/, where f eE+H, so that {dt + TBo)Ft = 0. Then 



Bo' 



Proof. To prove the first identity, note that 

= -i{BoN+ - N' Bo)-\Bod + d*Bo)Ft + i{BoN+ ~ N- Bo)-\Bod + d*Bo)Ft 

= t{BoN+ - N-Bo)-\{Bo - Bo)d + d\B^ - Bo))Ft 



since BqN~^—N Bq = BqN~^—N Bq. The vector part of this matrix is 
Furthermore 







Bo-Bo = 















0, 



which shows that 



{Bq — Bo)dFt = ttj^^^dF^'^ and d* {Bq — Bo)Ft = d*a^\j^Fl' 
using the mapping properties of d and d* from Remark 12. 6[ Thus 



-la 



dF^'^ + d*a^^^Ft 



)F}' 



since BoN~^ — N Bq = — a^x on normal vector fields. 

To prove the second identity, we multiply the Dirac equation {dt + TBo)Ft = by 
BoN+ - N-Bq to obtain 

(5oiV+ - N-Bo)dtFt - i{Bod + d*Bo)Ft = 0. 

The normal component of the vector part on the left hand side is 

—a^^dtFl'^ — i{a^\^dFl'^ + c^*('^iix-^/'^ + ^tiiii-^/'")) = 0. 

Here we have used the expression for the vector part of BqN^ — N^Bq above for 
the first term. For the other terms we recall from Remark 12.61 that the vector part 
of dFf is dFl'^ which is tangential, and that the normal vector part of (TBqFi is 
c?*(-Bo-^t)^'"- Therefore, after multiplying the equation with o^q, we obtain 

«x|i)^^/'^ - ^c?^o'^a-^a||^)F/'^ = dtFl'^ + id%^Fl'\ 

which proves the lemma. 

As in Lemma 13.11 we note that since I3q is a block matrix 

T^^ = = F + Bq^T*Bo, where F := -iNd. 

To prove Theorem 14.121 we shall need the following corollary to Theorem 14. 8[ 



□ 
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Corollary 4.15. Let i?o be the block matrix defined above, assume that WQtWojj < C 
and let v G Loo(R"; A) be a function with norm ||f ||oo < C . Then for all f G E^^Ti 
we have estimates 



(4.17) 
(4.18) 



\\\Q,£eTTlPt^\Fl^\)\\\ < (|||g,|||„, + l)||^||oo 
|||Q,fP^(F/'%)|||<(|||g,|||„, + l)||^|U 



where Ft := e ^^^'^ol f and Ff ' ={Ft,eo). The corresponding estimates for f G E ^^H 
also holds. 

We defer the proof until the end of this section, and turn to a lemma in preparation 
for the proof of Theorem 14.121 

Lemma 4.16. If f e E+ H and Ft := e'^^^'^olf, th 



en 



(4.19) 
(4.20) 



s . ^ ^ .ds 

-(.a,F,)- 



< 



\\\tdFt\\\ < 



^tFtlW < 



The corresponding estimates for f G E^^Ti also holds. 

Proof. The proof of (14.191) uses Schur estimates. Applying Cauchy-Schwarz inequal- 
ity, we estimate the square of the left hand side by 



dt 

— < 

t - 



* s ds 
t~ 



s ,, ^ ^ ,,nds\ dt 
t" .s t 



2 < 



where ip{z) := ze~'^L 

To prove fl4.20p . we note that iMbqTbq = d + d*B^ and write P^q, P^q for the two 
Hodge projections corresponding to the splitting 7i = N(rf) © N(d^^J. Note that 
these projections are bounded by Lemma [!^.21( i). We get 

tdFt = i-p\^MB,{tTB,Ft) = iP\,MBMTB,)f, 
from which fl4.20p follows. □ 

Proof that Corollary 4.15\ implies the estimate (4.14). By inserting / = Pj^" + (/ — 
PF") we write the left hand side in (14.141) as 

QtSPt^'iFt - /) + QtS{I - Pt^')Ft - QtS{I - Pt^°)f =: X,+X,- X3. 

For X3 we get from Theorem 14.81 the estimate |||X3||| < (|||Qi|||op + l)||^^||c 



smce 



/ - pfo = t2pp^^pBo ^ t^Y*^TPt^°. For the term X2 we use the first identity in 
Lemma [4. 141 to obtain 

X2 = QtSQfHT^^Ft = -QtSQf°{tdtFt) - iQtSQf'a^o'a^,{tdFl'^) 

- iQtS{QfHd%J{Fl'\a-la,^eo)) =: -X4 - zX, - tX^. 

We have the estimate |||X4||| < ||i^||oo|||tc^t-^t ||| ^ ll^l|oo||/||- For X5, we see from 
Remark 12.61 that dFl'^ is the vector part of dFt. Thus |||X5||| < ||£^||oo|||^^-Ft||| ^ 
||^||oo||/|| by ^M). To handle the term Xq we note that QfHd) 



-Bo 
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using Remark 13. 2[ Thus we obtain from Corollary 14.151 with v 



a||/a||xeo, the estimate |||X6||| < (|||(5t|||op + 

It remains to estimate the term Xi. To handle this, we separate the normal vector 
part as 

Xi = gi^P,^«(F/'%) - gi^P,^°(/i'°eo) + - g) =: Xj-Xs + X,, 

where Gt := Ft — P/'^ and g = f — /^'^. From Corollary I4.15[ with v = Cq, 
we get the estimate HlXylH < (|||(5t|||op + l)ll^||oo||/||- For the term Xg, we write 
Pf° = / - t^TT% - t^r*. FP^". From Theorem WE we obtain the estimate 

^ ^ -Bo -Bo ^ 

lll^slll ^ (IIIQtlllop + l)ll'^l|oo||/||- For the term Xg, we integrate by parts to obtain 
Xg = QtEPt'itdtGt) - QtEPt' [J^ sd^GJs^ =: Xio - Xn- 



We have the estimate |||Xio||| < ||£^||oo|||ic^t-^t||| ^ ||^^||oo||/||- For the term Xn we 
apply dt to the last expression for {dt + Tj^^)Ft in Lemma 14.141 and get 

(4.21) d'.Gs + d.T^F, = id^dsFl^K 

Thus 

Xn = -QtS{PSW (^J^ jisdsF,)^'^ 

+ iQMPSdl,) (£^(s5.P,^'ii)^^ =: -Xi2 + ^Xi3. 

Both IIIX12III and |||Xi3||| can now be estimated with ||£^||oo||/|| by f l4.19p since 
Pf"tT^^ = Qf° and Pf°td*^^^ are uniformly bounded by Corollary 12.231 This proves 
the estimate (Oil) . □ 

Proof that Corollary 4-15\ implies the estimate 1^4. 15 ). We write the left hand side 
in (14.151) . using integration by parts, as 

Qt£d I FJs = QtStdFt - QtSd [ sdsFJs =: Xi - X2. 



For Xi we have \\\QtS{tdFt)\\\ < \\S\U\f\\ by ^Ml- For X2, we write / = Pf'^ + 
(I -Pt^") and get 

X2 = Qt£itdPf') / -{sdsFs)— + QtStdQf' / sa,T^^P,ds =: X3 + X4. 
Jo 't ^ Jo 

Using that \\tdPf''\\ < C by Corollary E^S and fICTD shows that IIIX3III < ||f ||, 
To handle X4, we use the identity (I4.2ip . which gives 

X4 = QtStdQf" (- [ sd'^Gsds + i [ sd% 9,P,^'"rfs^ =: -X5 + iXg. 



For Xe we note that dQ^'^d*^^ = NdQ^^Nd^^ = -TQf'T*^^ = -T^Qf" = using 
Remark 13.21 and thus Xg = 0. To handle X5, we rewrite this with an integration by 
parts as 

= QtS{tdQf°){tdtGt) - QtStdQt^Gt + QMtdQf^^g) =: Xg - X7 + Xg, 
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where g = f - /^'^. Using Lemma EE] and that \\tdQf°\\ < C by Corollary E^Sl 
we get IllXelll < ||^||oo||/||. For Xg, we note that trfQf" = iNt'^TT*^^Pf'> . Thus we 

can apply Theorem iSl to obtain HlXglH < (|||Qt|||op + l)||^||oo||/||- 
We now write Gt = Ft — F^'^eo and get 

Xr = QtStdP^^tT^^Ft) - Qt£tdQf\Fl^\,) =: X, - X,^. 

Again noting that tdQf° = iNt'^TT*^ Pf° we obtain from Corollary 14. 151 with v = eo, 

the estimate |||Xio||| < (|||Qt|||op + l)ll^l|oo||/||- The term Xg remains, on which we 
use the first identity in Lemma [4.141 to obtain 



Xg = -QtStdPf^itdtFt) - iQtS{tdPf')aoo'a^^^{tdFt 



l,-L^ 



- tQt£itdPSd%,)iFl'\a-la,^eo)) =: -Xn - tX^^ - iX 



13- 



Using that prfPf"!! < C by Corollary ESHl shows that |||Xii||| < ||^||oo||/||. For 
Xi2 we see from Remark 12.61 that dP}'^ is the vector part of dPt. Thus |||Xi2||| < 
ll^lloolllt^i^tlll ^ ll^lloo||/|| by ^UM . To handle the final term X13 we note that 
dPP°d\ = dNPf'^Ndl = -FR^T*. = -rr*~ Pf° using Remark O Thus we 

obtain from Corollary 14.151 with v = a^^l^a^\j_eQ, the estimate |||Xi3||| < (|||Qf|||op + 
l)||^||oo||/||. This proves the estimate fHlHD . □ 

Proof of Corollary 4-15 Let / G E^J-L and consider the functions Pt = e~*'"^^ol/ and 
Gt ■= Pf»f. We note that P^ - Gt = iptiTB,)f, where ip{z) = e-\'\ - (1 + z"^)-^ e 
^{Sl). Thus it suffices to prove the estimate |||e^(G't'°)||| < (|||Qt|||op + l)||^||oo||/||, 
i = 1,2, for the two families of operators 

(4.22) el := Qt£t'TT*^P^'>M„ 

(4.23) e? := QtSPt^'M,, 

where My denotes the multiplication operator My{f) := vf. To this end, we write 

QliG]'') = (ej - 7M^)G'^° + llAG]^', 

where 7t(x) is the principal part of the operator family as in Definition 14.21 Using 
the principal part approximation Lemma 14.41 and Carleson's Lemma 14.61 we obtain 
the estimate 

\\mGl'')\\\ = ||ei||ofr|||tV(G,^'°)||| + H\\c\\NMtG't')l 

By Proposition 12.251 and Lemma [2.261 we have ||0t||off ^ ll^lloo- Furthermore, by 
Lemma [4.71 and Theorem 14.81 we have 

Il7lllc<ll|eiii|op + i|eiiioff<(|||g*|||op + i)||^||oo. 

For 92 we have used that 6^ = QtS{I - t^VV*^^P^° - t^F^^ FPf °)M^. 
To bound |||tV(G'J'°) |||, we note that the vector part of dGt is 

dG]'^ = zmci(eoGt'°) = -id{Gl'^) = -zV(G^°). 
Thus, similar to the proof of fl4.20p . it follows that 



m{G't')\\\<\mGt\\\<\\\tTs,Pt''\ 



< 
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Finally, to bound \\N4AtG]'°)\\ we write Gt = \{Ht + //-t), where Ht := (1 + 
itTBo)~^f- We now observe that the divergence form equation and estimate for hI'^ 
in Lemma [2.571 in fact holds for all / G ?i, by inspection of the proof. We get 

N,{AtHl'°){x) = sup / Hl'\z)dz 

\y-x\<t \jQ{y,t) 

< sup (J \Hl''iz)\''dzY' < (M(|/n(a:))i/^ 



t>0 \/B{x,rot) 

where Q{y,t) denotes the dyadic cube Q & At which contains y. Using the bound- 
edness of the Hardy-Littlewood maximal function on L2/p(R'^), we obtain 

\mAHi^')h < ii(M(i/n)^/ib = iiM(i/nii^;j < wmff; = 

A similar argument shows that ||iV,(A/fif ) || < ||/||, and thus ||iV,(AGj'°)|| < 

||/||. We have proved that |||ej(G,''°) ||| < (|||Qt|||op + l)||^||oo||/||, and therefore 
Corollary ma □ 

4.3. Proof of main theorems. We are now in position to prove Theorems 
and 11.41 stated in the introduction. 

Proof of Theorem \1.4\ Given a perturbation B'^ G Loo(R"; /^(a'^)) of the unper- 
turbed coefficients G L^(R"; C{a'')), we introduce B := I ® ■ ■ ■ ® B^ ® . . . ® I e 
Loo(R";/:(A)) and 5o :=/©... ©5^ ©...© / G L^iK""; C{A)) acting in all H. By 
Theorem l3.3( i). Tb^ satisfies quadratic estimates and by Lemma l^?T] and Theorem l4.8l 
there exists e > such that we have quadratic estimates 

IIIQf /III ^ 11/11, whenever \\B - 5o||oo <e, feH. 

Therefore, by Proposition 12.321 and Proposition 12.421 we have when ||i? — i?o||oo < 
e/2 well defined and bounded operators Eb = sgn(TB) which depend Lipschitz 
continuously on B, i.e. 

\\Eb,-EbA\ < C\\B2 - BiW^, when ||S. - Soil < e/2, i = l,2. 

To prove that (Tr- i?^a;^) is well posed, note that by Lemma [2.511 it suffices to show 
that X — EbNb is invertible since then in particular X — EBkNBk = (X — EBNB)\'^jk is 
invertible. Here the spectral parameter is A = (a + l)/(a — 1) and a := /a~ . By 
Theorem l3.3( ii). the unperturbed operator X — Eb^Nb^ is invertible when A^ + 1 7^ 0. 
For the perturbed operator we write 

A - EbNb = (A - Eb,Nbo) (/ + pipj(A - NB^iEB^NBo - EbNb)) . 

Here WEb^Nbo — EbNb\\ < ||i? — i?o||oo since we clearly have Lipschitz continuity 
\\Nb2-NbA ^ IIS2-S1II00. It follows that A-EbA^b is invertible when ||S-5o||oo < 
C|A^ + 1|. Under the assumption ||i? — i?o||oo < e/2 with e small we can replace 
A2 + 1 = 2(^2 ^ ly^^ _ 1^2 ^-^^ ^2 ^ ^_ 

This proves that for each boundary function g there exists a unique solution / = 
f^ + f~ satisfying the jump conditions in (Tr-i?^a;^) and ||/^|| + ||/~|| ~ ||/|| ~ ||5'||- 
The norm estimates for in Theorem 11.41 now follows from Lemma I2.49[ Finally 
it follows from Proposition 12.421 that the solution operator 

F^{t,x) = 2e^*l^^lE|((a+ + a-)EB - {a+ - a-)NBy^g{x) 
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depends Lipschitz continuously on B. Indeed, from above it suffices to prove Lips- 
chitz continuity for e^*'^^'. For the trippel bar norm, we here take 4'{z) = ze^^^^ in 
Proposition I2.42[ and for ||-Ff'^||2 we take h{z) = e^*'^'. This completes the proof of 
Theorem O □ 

Before turning to the proofs of Theorems 11.11 and 11.31 we note some corollaries of 
Theorem 11.41 First, if we let = 1 in Theorem 11.41 then it proves that the following 
Neumann-regularity transmission problem is well posed for small Lqo perturbations 
A of a block matrix Aq. 

THE TRANSMISSION PROBLEM (Tr-Aa±). 

Let G C be given jump parameters. Given scalar functions : R" C 
with V^i) e L2(R"; C") and G L2(R"; C), find gradient vector fields F^{t,x) = 
Vt,xU^{t,x) in R!^+^ such that e C\R±]L2{W; C"+^)) and F± satisfies fOD 
for ±t > 0, and furthermore lim^^ioo = and limt^o± = in -^2 norm, 
where the traces satisfy the jump conditions 



a-V.[/+(0,x) -a+V,.t/-(0,x) = V,.7/^(x), 
a+f^(0,x)-«-f^(0,x) = <Pix), 

where \/xU^{0,x) = f^{x) and = {Af^,eo) denotes the conormal derivative. 

Secondly, Theorem 11.41 give perturbation results for the following boundary value 
problems for /c-vector fields. 

THE NORMAL BVP (Nor-B^). 

Given a fc- vector field g G Ti^, find a A;- vector field F{t, x) in R^+^ such that Ft G 
C^(R+; L2(R"; A^)) and F satisfies (Ol) for t > 0, and furthermore limi_oo Ft = 
and lim^^o Ft = f in L2 norm, where / satisfies 

eo J (EV) = eo J (S^^) onR" = 9R!;:+\ 

THE TANGENTIAL BVP {Tan-B''). 

Given a fc- vector field g G find a /e- vector field F{t,x) in R^^+^ such that Ft G 
C^(R+; L2(R"; a'')) and F satisfies ( Ol) for t > 0, and furthermore lim^^oo Ft = 
and lim^^o Ft = f in L2 norm, where / satisfies 

eoA/ = eoA^ on R'^ = (9R!^+\ 

Corollary 4.17. Let Bq = Bq{x) G Loo(R"; '^^(a'^)) be accretive and assume that Bq 
is a block matrix. Then there exists e > depending only on the constants ||i?Q ||oo and 
K^k and dimension n, such that if B^ G Loo(R'^; >C(A^)) satisfies \\B^ — -BqIIoo < ^, 
then the normal and tangential boundary value problems (Nor-B^) and (Tan-B^) 
above are well posed. 

Proof, (i) For (Nor-S'=) in R++\ we let a" = and a+ = 1 in {Ti-B^a^). Then we 
obtain two decoupled jump conditions 

-eoA/~ = eoA^f 
e,^{B^f+) = eoAB'g). 

Discarding the solution F~ , we obtain a unique solution F = to (Nor-i?^). 



FUNCTIONAL CALCULUS AND COMPLEX PERTURBATIONS 



63 



(ii) For {Tan-B'') in R++\ we let a' = 1 and a+ = in {Ti-B^a^). Then we 
obtain two decoupled jump conditions 

/ eoA/+ = eoA^f 

Discarding the solution F~, we obtain a unique solution F = to (Tan-i?*^). □ 

Note that well posedness of (Neu-A) and (Reg- A) for Theorem ll.l( b) is the special 
case = 1 of Corollary 14.171 as well as the special cases (a"*", a~) = (1, 0) and (0, 1) 
respectively of {Tr- Aa"^). 

Proof of Theorem M.^ Let Aq be such that Tb^ has quadratic estimates in ?i, where 
Bq = I®Ao®I®...®L Thus by Lemma lO and Theorem HE] there exists e > 
such that we have quadratic estimates 

IIIQf /III ^ 11/11, whenever \\B - B,\U <e, feH. 

Therefore, by Proposition 12.321 we have when ||i? — i?o||oo < e well defined and 
bounded operators Eb = sgn(TB). With Lemma [2.46^ these restricts to bounded 
operators Ej^ in TC^. In particular f Eli} can be decomposed as / = /"*" + /", where 
f^ := E^f and ||/|| ^ ||/+|| + ||/"||. Moreover, Lemma and Proposition 12.561 
proves the stated norm equivalences for F^. □ 

Proof of Theorem \l.l\ Given a perturbation A G Loo(R"; C{/\^)) of the unperturbed 
coefficients Aq G Loo(R"; /^(A^)), which we assume are either of block form, real 
symmetric or constant, we introduce B := I®A®I®...®I G Loo(R"; 'C.{A)) 
and Bo := I ® Ao® I ® . . . ® I e Loo(R"; C{A)) acting in all H. That Tbo satisfies 
quadratic estimates follows from Theorem 13. 3( i). Theorem 13.121 and Proposition 13.51 
respectively. Theorem 11.31 now shows that we have quadratic estimates for Tg when 
ll-B — -BqIIoo < In case Aq is a block matrix, we note that this result follows 
already from Theorem 14. 8[ By Proposition 12.321 and Proposition 12.421 we have when 
\\B — -BqIIoo < e/2 well defined and bounded operators Eb = sgn(TB) which depend 
Lipschitz continuously on B, so that 

\\Eb, - EbA < C\\B2 - B^W^, when \\Bi - BoW < e/2, t = 1,2. 

To prove that (Neu-A) and (Reg-A) are well posed, note that by Lemma 12.531 it 
suffices to show that / ± E^Na : —>■ are invertible, where E^ = EbI-^i and 
Na = NbI-^i- By Theorem I3.3( ii). Theorem 13.121 and Theorem 13.61 respectively, the 
unperturbed operators I ± EaqNaq '■ n Ti.^ are invertible. For the perturbed 
operator we write 

I ± EaNa = (/ ± EaM (/ ± (/ ± EaM-\EaNa - EaM) • 

Here H-EyiA^yi — EaqNaqW < ||A — Ao||oo since we clearly have Lipschitz continuity 
\\Na2 — ^Aill ^ 11^2 — ^i||oo- It follows that / ± EaNa are invertible when ||y4 — 
^lloo < e'- 

The well posedness of (Neu-'"-y4) is a consequence of Proposition 12.521 since our 
hypothesis is stable when taking adjoints A ^ A* . Alternatively, we can replace 
Naq with above, proving that I ±EaN is an isomorphism, using Theorem 13. 3( ii). 
Remark 13.131 and 13.71 This proves that for (Neu-A), (Reg-A) and (Neu-'"-y4) and 
each boundary function g (being (p and Vxip respectively), there exists a unique 
solution / = E^Tv satisfying the boundary condition and ||/|| ~ \\g\\. Lemma [2.491 
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and Proposition 12.561 proves that the stated norms of / and F{t,x) := (e-*l^^l/)(a;) 
are equivalent. 

The well posedness of (Dir-A), as well as the first three norm estimates, follows 
from Lemma [2.551 and (Neu-'"-yl). To show that |||tVa;?7t||| ~ we consider 

the gradient vector field Gt = Vt,xUt as in the proof of the lemma. From (Reg-A) 
and (Neu^-v4), it follows that for all t > 0, we have \\dtUt\\ = ||A^"Gt|| ^ ll^^+GiH = 
llVa^^f/J, from which the square function estimate for VxUt follows. To show that 
||n|| ^ ||A^=K(f/)||, we consider the vector field Ft of conjugate functions from the 
proof of the lemma. Proposition 12.561 shows that ~ ||/|| ~ ||A^*(F)|| > ||A^=K(f/)||. 
Moreover, the proof of the reverse estimate < ||A^*(f/)|| is similar to the proof 

of 11/11 ~ ll^*(-^)ll i'^ Proposition I2.56[ using the uniform boundedness of Vt- 

Finally we note that the solution operators for (Neu-A), (Neu-'"-A), (Reg- A) and 
(Dir-y4) are 

Ft = 2e-'\^^\EA - iV^)-i(aooVeo), = 2e-*l^^l(E^ - iV)-i(0eo), 

Ft = 2e-*l^^l(E^ + iV)-i(V.V'), Ut = 2(e-*l^-^l(E^ - NY\ue^), eo). 

The Lipschitz continuity of m i— > t/ is a consequence of the corresponding result 
for (Neu-'"-A). For the norms sup^^g \\Ft\\ and |||t9i-Ft|||, Lipschitz continuity for the 
solution operators follows from Proposition 12.421 as in the proof of Theorem II. 4[ It 
remains to show Lipschitz continuity for the norm H-FHa" = ||^*(-^)||2- To this end, 
we consider 

f{x) ^ Fa. = {e-'\^'^\ElJ){x) : ^ X 

and the truncations f{x) F^^(t,x) = Xkit)FA.{t,x), where Xk denotes the char- 
acteristic function for {l/k,k) as in the proof of Lemma l2.41( iii). We claim that 
it suffices to show that, for each fixed k, the operator f{x) ^ F^ {t,x) : Ti.^ X 
depends holomorphically on z. Indeed, using Schwarz' lemma as in Proposition 12.421 
we obtain the Lipschitz estimate 

\\Fl,- FX^\U <C\\A,- A,\U\f\\„ 

uniformly for all /c, since HF'^H^f < < C\\f\\ by Proposition l2.56[ Furthermore, 

by the monotone convergence theorem we have 

\\FX- F\^\U / \\Fa,- FaAx. k^oo, 

so the desired Lipschitz continuity follows after taking limits. 
To prove that f{x) ^ F\ {t,x) is holomorphic, we note that 



JR" \t>0 J\s-t\<cot J\y-x\<cis/(l-co) ^ J 

< f it I \F\s^y)\-^-^\dx^ I WF^Wlds^ 



JK" \Jl/k J\y-x\<cis/{l~co) ^ / Jl/k 

for fixed k, where in the last step we use that 1/k < < k. Since f{x) t— > F^{t, x) : 
L2(R" X [l/k,k)) is holomorphic by Lemma l2.41( ii) and the embedding 
L2(R" X {1/k, k)) "—>■ X is continuous and independent of z, it follows that f{x) i— >■ 
F^{t,x) : Ti.^ — > A* is holomorphic for each fixed k. This completes the proof of 
Theorem 11.11 □ 
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